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ABSTRACT 


The subject of this investigation is to determine the transformation 
parameters (three rotations, three translations and a scale factor) between 
two Cartesian coordinate systems from sets of coordinates given in both 
systems. The objective is the determination of well separated transformation 
parameters with reduced correlations between each other, a problem especially 
relevant when the sets of coordinates are not well distributed. The above 
objective is achieved by preliminarily determining the three rotational para- 
meters and the scale factor from the respective direction cosines and chord 
distances (these being independent of the translation parameters) between the 
common points, and then computing all the seven parameters from a solution 
in which the rotations and the scale factor are entered as weighted constraints 
according to their variances and covariances obtained in the preliminary solu- 
tions. 

Numerical tests involving two geodetic reference systems were performed 
to evaluate the effectiveness of this approach as follows: 

(a) A non-constrained solution for general transformation for the seven 
parameters (including the three translations and scale factor). 

(b) A constrained solution for general transformation for the seven 
parameters utilizing the three rotations with their statistics as 
constraints. 

(c) A constrained solution for general transformation for the seven 
parameters using the three rotations and scale factor with their 
statistics as constraints. 

The above schemes were then separately repeated for each of the following 
three cases: 
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(i) Using the full variance -covariance matrix between coordinates 
of the geodetic reference systems. 

(ii) Using only a (3 X 3) banded diagonal variance-covariance matrix, 
thus assuming no correlation between coordinates of any two 
points within the system. 

(iii) Using only variances for the coordinates, thereby further omitting 
the correlation between the three coordinates of any one point in 
the system. 

In the case of seven parameter general transformation, the best estimates 
were obtained using full variance-covariance matrix and constraining three 
rotations and the scale factor, case (c) and (iii) above. The improvement in 
correlation between translations and rotations was more significant compared 
to between translation and scale factor. 
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1. INTRODUCTION 


During the last twenty-five years with the availability of computer technology 
and its phenomenal growth in basic hardware and core storage capacity and the 
exceptional increase in a computer’s ability of solving problems in lesser and 
lesser time, a trend has set in to analyze the problems in geodesy and photogram- 
metry more and more in three dimensional space rather than to follow traditional 
concepts. 

Further, the advent of artificial satellites and their subsequent use in geodesy 
made it possible to obtain Cartesian coordinates of points on earth surface. 

Several projects involving satellite-networks of continental or global extent 
were begun and at present they are in varying stages of completion. Many new 
solutions have recently come out, each delineating its own reference system. These 
systems in reality should differ from each other only in having different origins, 
sets of axes or scale. 

Thus, the relationship between any two such reference systems (e. g. , UVW 
and XYZ) would generally consist of seven parameters — three translations (AX, AY,AZ) 
between the two origins, three rotations (ci;0,f) of the Euler's angle type between the 
two sets of axes and the scale factor (As), if any (Figure 1). 
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The mathematical model to be used in the computations of the above seven 
parameters from a least squares solution may be written in the following form 
fBarlokas, 1969; Bursa, 1965; Wolf, 19631: 
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where "i" denotes any point common to both the systems. The three angles oo, 
0,and € of the Euler type correspond to small rotations about the Z, Y and X axes 
respectively — the positive direction of rotations taken in counter clockwise mode, 
when viewed from the end of the repsective axes towards the origin. It may be 
worth while to mention here that the station coordinates in both the systems 
(Ui , Vi , ;and , Z^) are treated as observations in the above model. 

The above equation written in matrix notation can then be modified into the 
observation equation below [Uotila, 1967]: 


BV + AX + W = 0, 
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while V and X represent the residuals to the observations and corrections 
to the parameter estimates, respectively. Hence, collecting all the 
matrices as above, pointwise in the systems, the observation equation 
becomes: 
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(3) 


Defining the geodetic reference systems on the assumption that the 
Laplace -condition has been enforced throughout the network (which implies 
that the axes of the reference ellipsoid are parallel to the conventional 
earth-fixed axes), many experiments have been made in recent times 
to determine the seven transformation parameters in relating the different 
geodetic systems to each other using an observation equation of type (3) 
TLambeck, 1971; Marsh et. al., 1971]. 



However, in the above general transformation, if the geodetic reference 
systems are properly oriented through the Lapl ace-condition, the three rotations 
arising due to the improper relative orientation of the systems are generally never 
more than a few seconds of arc, while translations may amount up to 200 to 300 
meters. Also, due to the presence of high correlations between the rotations, 
the scale factor and the translations, satisfactory independent estimates for these 
parameters are difficult to obtain from a combined general solution using equation (3). 

This investigation separates the determinations of the rotations and the scale factor 
(from that of the translations) for subsequent use as constraints in a combined general 
solution. 


2. THE INDEPENDENT DETERMINATIONS OF 
ROTATIONAL AND SCALAR PARAMETERS 

2. 1 Determination of Rotations 
2. 1. 1 Mathematical Model 

The mathematical model used in this study is as follows [Bursa, 1966]: 

- t 4^^ + C0 + 0sinTjk ^ tanS/ic^^ - CcosT^^ tanSii’'^ 

+0cosT/k^ +esinTik^ 

where Tm and 6uc are defined as the geodetic hour angle and declination of the 
(i direction of the observed point at station and the observer at i'^*‘ station. 

The indexes (1) and (2) denote the two systems with the transformation proceeding 
from system #1 to system #2. 

If Aiic , Bik , Ciic are taken to denote the direction cosines of the (i -k)^'’ line of 
length Ri^, then for the first (UVW) system one gets: 
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In the above relations (4) through (6) the elements of translation do not enter 
the picture A similar set of relations as per (5) and (6) can be established for the 
second (XYZ) system. 


2.1.2 Observation Equations 

The mathematical model (4) then, for each (i-k)^" line, yields the following 
generalized form of observation equations fUotila, 1967]: 
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Using the conventional weight matrix P for the coordinates of points included in 
the transformation (see section 2. 1. 3), and the principle of least squares by making 
V PV as minimum, the equation (7) is then solved for correction vector (co,ip, c) and 
for the variance-covariance matrix (Eoo0c) of the three parameters. 
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2.1.3 Weights 

Using the variance -covariance matrices EX and EU in respect of i*'*' and k'’*' 
points for the XYZ and UVW systems, the variance-covariance matrices Et 6 
the two systems of coordinates can be computed through propagation of errors 
[Uotila, 1967]. 

Two distinct cases would arise here. Firstly, when in addition to correlation 
between X, Y, Z-coordinates of any point, the correlation between the coordinates 
of one point to others is also considered. In such a case, the necessary relation will 
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Secondly, ignoring the correlations between the coordinates of different 
points within a system, equation (8) can be modified as under: 
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In the equations (8) and (9), LUi and LUic correspond to i^** and point of the 
first system and can be either full (3 x 3) matrices with covariances between the 
three coordinates of a point, or may contain variances for U, V and W in a diagonal 
form only. However, in the case of covariances (LUuc) between the points being 
included, the matrix in equation (8) would be a full (6 X 6). 

Obtaining similarly , the combined variance-covariance matrix, to 
be used with equation (7), is given by: 
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It may be noted here that the matrix P is always in 2 x 2 banded diagonal 
form. 


2. 2 Determination of Scale Factor 
2.2.1 Mathematical Model 

The scale factor between the systems #1 and #2 would be given as follows: 

p (2 ) 

^Suc = “ 1 (11) 

“Ik 
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where 
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2.2.2. Weights 

Using the variance-covariances matrices EX and EU for the coordinates of 
ij’’ and points in the two systems included in the transformation (section 2, 1. 3), 
a variance a®g is established for the scale factor through error propagation. Two 
cases similar to equations (8) and (9) would arise according to the case when full 
variance-covariance matrix between different points within the system is considered 
or not. 

The matrix G for the scale factor determination is 
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Hence, 
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where the full (12 x 12) matrix would become a (3 X 3) banded diagonal matrix in 
case EUik and EXi^ are zero, 1. e. , covariances are not considered. The 
complete (12 x 12) matrix would assume a diagonal pattern when only variances 
are used for station coordinates. 

Using the value of Asjjt and from equations (11) and (12), the value 
for weighted mean and its variance for the transformation under investigation 
is established as given below [Hirvonen, 1971]: 
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where 
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Ik 
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l/cTa,i^ and fwu,] denotes the sum of all such weights. 
Total number of scale factor values used in the sample. 


3. BRIEF DISCUSSION ON THE FORTRAN PROGRAM 

Appendix I gives the complete computer program for obtaining the constrained 
or non-constrained solution fOr seven parameters. With appropriate coding 
non-constrained solutions for three parameters (AX, AY and AZ) and scale 
factor ^ can also be obtained. 

The input coordinates can either be Cartesian or geodetic (ellipsoidal) with 
35 as the maximum number of points in each system. However, the matrices 
can easily be re-dimensioned to accomodate more points when required. The 
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program is self-explanatory with regard to definition of various option codes for 
input, type of solution and inclusion of correlation data, etc. 

The broad basic divisions of the program are as under: 

(a) Main Program ; This section takes as input the various options in 
input/solutions, coordinates of points, rectangular or ellipsoidal, and 
semimajor axis and flattening of the ellipsoid used, if required. It 
then prints out the two sets of coordinates used for checking purposes. 

The various options of input/solutions have been designated in the 
program as KCODE e.g. , KCODE (1) refers to number of common 
points involved in the transformation. A complete list with necessary 
explanatory remarks has been included in the beginning of the program. 

(b) Subroutine "EULERS” ; This subroutine first reads the variance-covariance 
matrices of the station coordinates, with or without correlation, and 

then sets up matrices A, W and P to be used for the solutions of three 
rotations through direction cosines (equation (7)). 

The subroutine writes up the variance-covariance matrices for 
the coordinates on the disk and stores the estimates for co,j/) and c, 
and their variance-covariance matrix [Lcoi/)€] in the common block for 
subsequent use. 

(c) Subroutine "SCALE" ; This subroutine computes the weighted mean 
value for scale factor ^s and its variance by direct chord comparison 
independent of other transformation parameters (equations (13) and (14)), 

(d) Subroutine "TFORM" ; This subroutine solves for a general transformation 
(equation (3)), utilizing the common block core memory for coordinates 

of points and variance-covariance matrices from the disk. 

The matrix M'^ to be utilized for generating normal equations is 
computed by calling another subroutine "SETUP", 

NOTE ; In case the solution is required ONLY for three translation 
or three translations and scale factor, KCODE (3) is coded as "0" 
and then subroutine "EULERS" is skipped by the program. 
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(e) Subroutine "CSTRNT" : This subroutine uses the results of subroutines 

SCALE and EULERS as constraints with their appropriate statistics 
and computes for a constrained solution of seven parameters. The 
results are returned to subroutine TFORM for printout. KCODE (11) 
refers to the option whether 3 or 4 parameters are to be constrained. 

(f) Subroutine "RESIDU" ; This subroutine computes the residuals vector 
V for observations i. e. , the station coordinates used in the program. 
The residuals are printed station wise for both systems #1 and #2. 

In the computer program, the storage mode used for major computation is in 
vector form for increased flexibility and saving of core storage. 

Appendix II gives a typical set of Job Control Cards (JCL). 
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4. NUMERICAL EXAMPLE 


The above transformation models were used to study the relationship between 
the transformation parameters and obtaining their best estimates by minimizing 
correlation for the following two reference systems: 

(i) System MPS-7, [Mueller and Whiting, 1972]. 

(ii) System NA-9, [Mueller et. al. , 1972], 

Using the same set of thirty common stations of the above two systems, the 
following solutions were obtained during the investigation; 
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©Note: The constraints for these solutions (rotations and/or scale factor) 

with their statistics were computed independently of the translation parameters 
(subroutine EULERS and SCALE of the Fortran IV program). 

Two solutions in full have been appended in the report as specimens in Tables 
1 and 2 as under: 

Table 1: Sample printout of the solution for three rotations (CjO,0 , €) and 

scale factor (As), using full variance-covariance matrix. 

Table 2: Sample printout of the constrained seven parameter general 

solution between NA-9 and MPS-7 with three rotations and 
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TABLE 1 


Sample Printout of the solutions for three rotations as 
parameters and the scale factor, using full variance- 
covariance matrix. 
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TABLE 1 


SOLUTION TUK »3» ROTATION PAKAMETbKS 
(FRUM DIRECTION COSINES — UNITS SECONDS OF ARC) 
(USING FULL VARIANCE-COVARIANCE MATRIX) 

OMEGA PSl tPSiLON 

C .io03791P+p0 -0.31>20iA5a--01 -0 . 2173oiuD+vO 

variance - COVARIANCE MATRIX 

MC2- 

0. 167?.3661D-0? 

0.A0622267D-03 
— C» 9i 7 677 6AD— 03 

0.iCOoOOOOD>Oi 
O.26277933D+0C 
-0.^393i5QiD+C0 


0*406232870-03 
0.123179910-02 
-0.48803740D-C3 
COEFFICIENT OF CORRELATION 


-L • 93 7 67 7 0 ‘I U “0 3 
-U.4880374UD-03 
C. 271919360-02 


0. 282779330+00 
0 , 100000000+01 
-0.266663210+00 


-0.439315010+00 
-0.266663210+00 
0. 1 00 00 GOOD +01 


SOLUTION FOR SCALE FACTOR 


(FROM CHORD COMPARISON) 


SCALE FACTOR 
(10.D45) 


VARIANCE 
(10.D411 ) 


5.10 


0.06 
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TABLE 2 


Sample printout of the constrained seven parameters 
general solution, using full variance-covariance matrix 
{case (c)/(iii)). 
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TABLE 2 


SCALE FACTOR ANU ROTATION PARAMETERS CONSTRAINED 


SOLUTION FUR 3 TRANSLATlQNy 1 SCALE AND 3 RGTAT1CB4 PARAMETERS 

♦ 

(USING FULL VARIANCE-COVARIANCE MATRIX) 

OX DY OZ OL OMEGA PSI EPSILON 

METERS METERS METERS 110 , 0 + 5 ) SECONDS SECONDS SECONDS 

-45.38 171.94 187.44 5.14 0.17 -0.04 -0.22 

VARIANCE - COVARIANCE MATRIX 


M02= 0.84 


0»I760^01 

0.2500+00 

0.4530+00 

-0.3I0D-07 

0.1260-06 

0. 7780-07 

-0.6520-07 

0*2S0D+00 

0.2280+01 -0.3220-01 

0.2430-06 

0.5310-07 

0.238D-07 

-0.1240-06 

0.453D+C0 

-0.3220-01 

0.2060+Ul 

-C. 1490-06 

0.6150-07 

0.2220-07 

-0.177D-06 

-0#3l0D-07 

0. 2430-06 -0. 1490-06 

C. 4410-13 

-0.3230-17 

-0.2980-16 

-0.1270-16 

0.126D-06 

O.55XD-07 

0.615D-07 

-0.3230-17 

0.2250-13 

0.5250-14 

-0.125D-13 

0.7780-07 

0.2380-07 

0.2220-07 

-0.29SD-16 

0.5250-14 

0.1670-13 

-0.6340-14 

-0.8S2D-07 

-0.1240-06 -0.177D-06 

-0.12 70-16 

-C. 1250-13 

-C.654D-14 

0.364D-13 


COEFFICIENTS OF CORRELATION 


0.1000+01 

0. 1250^00 

0.2380+0C 

-0.1110+00 

0.6350^00 

C.454U400 

-0.337D+C0 

0. 1250+00 

O.lOOD^Ol 

-0.1490-01 

0.7650+00 

0.24404^00 

C. 12204-00 

-0 .429 04-00 

O.238O+G0 

-0.1490-01 

O.lCOO+Ul 

-C. 4930+00 

0.286D4>00 

0.12CO'400 

—0 .64604-00 

-0. 1110+00 

C.7650^00 

-0. 4930+00 

0.1000+01 

-0.1030-03 

-O.llOD-02 

-0.3170-03 

0.6350+00 

0.244D+00 

0.2860+00 

-n.i03D-03 

O.ICOD-KM 

0.27104-00 

-0.436L^OO 

0.4540+00 

U.1220'^00 

0.1200+00 

-0.1100-02 

0.27104-00 

0.10004-01 

-0.263D4-C0 

-0.3370+00 

-0. 4290^00 

-0.6480+00 

-<». 3170-03 

-C. 43604-00 

-D.263D4-00 

O.lOOL'4-01 
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scale factor as constraints, using full variance-covariance matrix 
(case (c)/(iii)), 

A summary of the results for cases (a) through (c) and (i) through (iii) are 
presented in the following tables: 

TABLE 3 gives the results for three rotations, as obtained independently 
of translations and scale factor from direction cosines, for cases (i) through 
(iii). 

TABLE 4 gives the results for the scale factor, as obtained by direct chord 
comparisons independent of other transformation parameters, for cases (i) 
through (iii). 

TABLE 5 gives the results for the constrained and non-constrained seven 
parameters general transformation solutions (cases (a) through (c) and (i) 
through (iii)). 

TABLE 6 gives the comparative study of the results for seven parameters 
general transformation solutions as regards correlation between translations 
and rotations/scale factor, using different variance-covariance matrices (cases 
(i) through (iii)). 

TABLE 7 gives the comparative study of the results for seven parameters 
general transformation solutions as regards correlation between translations 
and rotations/scale factor, using different constraints (cases (a) through (c)). 
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TABLE 3 


Three Rotation Parameters from Direction Cosines 


NA-9~MPS-7 



Using Variances 
Only 

Using (3X3) Banded 
Diagonal Variance- 
Covariance Matrix 

Using full Variance- 
Covariance Matrix 

Case 

m 

^ (ii) 

m 


0. 17 ± 0. 05 

0. 17 ± 0. 04 

0. 17 ± 0. 04 

0 (") 

0. 04 ± 0. 04 

-0.02 ±0.04 

-0.04 ±0.04 

M") 

-0, 20 ± 0. 06 

-0. 24 ± 0. 05 

-0. 22 ± 0. 05 

CTo 

1. 15 

1.30 

1.36 


TABLE 4 

Scale Factor From Chord Comparison 
NA-9-MPS-7 



Using Variances 
Only 

Using (3X3) Banded 
Diagonal Variance- 
Covariance Matrix 

Using full Variance- 
Covariance Matrix 

Case 

(i) 

(ii) 

(iii) 

.As(X10®). . 

5.46 ± 0. 24 - 

5. 37 ± 0. 24 

5. 18 ± 0.24 
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TABLE 6 


Comparative Study of Correlation Coefficients 
Between Transformation Parameters 
(Using Different Variance-Covariance Matrices) 


Case (i): USING VARIANCES ONLY 



Non-C onstrained 

Constrained Solutions 


Solution 

3 Rotations 

3 Rotations and 
Scale Factor 

Case 


(a) 


(b) 

(c) 

'r'T’Tr-' — Translations 
Rotations — -LI 

and Scale. Factor 



AZ 

LX 

LY 

D 

LX 

LY 

LZ 

U) 

0.88 










0.63 





0.08 

0.51 



c 

-0.47 





Ba 

-0.40 


-0. 73 







■ 

-0. 10 


-0.44 


Case (ii): USING (3 X 3) BANDED DIAGONAL 
VARIANCE-COVARIANCE MATRIX 



Non-C onstrained 
Solution 

Constrained Solutions 

3 Rotations 

3 Rotations and 
Scale Factor 

Case 

(a) 

(b) 

(c) 

Translations 

Rotations ' 

and Scale Factor 

AX 

ay 

LZ 

AX 

AY 

D 

AX 

AY 

LZ 

CO 

0.83 


0. 33 

0.58 



0, 62 

0.24 

0. 27 


0.54 


0. 13 

0.38 



0.40 

0. 12 

0. 13 

€ 

-0.45 


0, 80 

-0.32 



-0. 34 

-0. 44 

-0. 66 

As 

-0. 15 


-0.56 

-0. 36 




0.76 

<?5 

O 

1 
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TABLE 6 (Continued) 


Case (iii): USING FULL VARIANCE -COVARIANCE MATRIX 



Non -Constrained 
Solution 

Constrained Solutions 

3 Rotations 

3 Rotations and 
Scale Factor 

Case 


(a) 

(b) 

(c) 

TrT’TT'— — _ Translations 
Rotati ons 

and Scale Factor ^ 

AX 

AY 

AZ 

AX 

AY 

AZ 

AX 

— 

AY . 

m 

CO 

0. 83 


0.33 

0. 60 


jm 

0. 64 

0.24 

0,29 


0.54 


0.13 

0. 43 



0. 45 

. 0.12 

0.12 

c 

-0.45 

-0.51 

-0.80 

-0.32 



-0.34 

-0.43 

-0.65 

As 

-0. 15 

0.84 

-0.56 

-0.36 


-0.89 

-0. 11 

0.76. 

-0.49 


- 21 - 
























TABLE 7 


Comparative Study of Correlation Coefficients 
Between Transformation Parameters 
(Using Different Constraints) 


Case (a); NON-CON STRAINED SOLUTION 



Using Variances 
Only 

Using (3X3) Banded 
Diagonal Variance- 
Covariance Matrix 

Using Full 
Variance- 
Covariance Matrix 

Case 

(i) 

(ii) 

(iii) 

Translations 

Rotations — 

and Seale Factor ' — — ~ 

AX 

AY 

Az 

AX 

AY 

D 

AX 

AY 

AZ 

00 

0.88 

0.40 

0.43 

0.83 

0. 27 

0.33 

0. 83 

0.27 

0.33 


0. 63 

0. 19 

0. 13 

0.54 

0. 11 

0. 13 

0.54 

0.11 

0.13 

c 

-0.47 

-0.67 

-0.88 

-0.45 

-0.51 

0. 80 

-0.45 

-0.51 

0.80 

As 

-0.10 

0.74 

-0.40 

-0. 15 

0.84 

-0.56 

-0, 15 

0.84 

-0.56 


Case (b): CONSTRAINED SOLUTIONS 
(CONSTRAINTS: 3 ROTATIONS) 



Using Variances 
Only 

Using (3X3) Banded 
Diagonal Variance- 
Covariance Matrix 

Using Full 
Variance- 
Covariance Matrix 

Case 

(i) 

(ii) 

(iii) 

Tr an si ati ons 

Rotations 
and Scale Factor 

AX 

AY 

AZ 

AX 

AY 

AZ 

AX 

AY 

m 

w 

0.68 

0.14 

0.22 

0.58 

0.09 

0.14 

0,60 

0.09 

0. 15 


0.49 

0. 07 

0.08 

0.38 

0. 04 

0.08 

0.43 

0.04 

0.07 

€ 

-0.38 

-0.23 

-0.45 

-0.32 

-0.16 

-0.34 

-0, 32 

-0, 16 

-0, 34 

As 

-0.29 

0.95 

-0.83 

-0.36 

0. 97 

-0. 89 

-0,36 

0.97 

-0.89 
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TABLE 7 (Continued) 


Case (c); CONSTBAINED SOLUTIONS 


(CONSTRAINTS: 3 ROTATIONS AND SCALE FACTOR) 



Using Variances 
Only 

Using (3X3) Banded 
Diagonal Variance- 
Covariance Matrix 

Using Full 
Variance- 
Covariance Matrix 

Case 

(i) 

(U) 

(iii) 

^r~T'rr- — Translations 
Rotations - 

and Scale Factor 



AZ 



AZ 

AX 

AY 

m 

CO 

0.71 

0.32 

0.35 

0.62 

0.24 

0.27 

0.64 

0.24 

0.29 

0 

0.51 

0. 14 

0. 13 

0. 40 

0.12 

0. 13 

0. 45 

0.12 

0.12 

c 

-0. 40 

-0.51 

-0. 73 

-0. 34 

-0.44 

-0. 66 

-0. 34 

-0.43 

-0. 65 

As 

-0. 10 

0.72 

-0.44 

-0. 11 

0.76 

-0.49 

-0, 11 

0. 76 

-0.49 






















5. CONCLUSIONS 


The comparison between different columns of Table 3 shows that the estimates 
for three rotation parameters remain more or less the same, but that their standard 
deviations show some improvement as we proceed from column 1 (variances only) 
to column 3 (full variance-covariance matrix). However, in the case of scale factor 
(Table 4) the estimates for As indicate a definite trend while standard deviation 
remains constant. 

In the case of seven parameters general transformation (Table 5) the comparisons 
among different columns indicate a definite overall improvement in all parameter 
estimates. The best estimates were obtained in the solution using full variance- 
covariance matrix and three rotations (u),0,c) and scale factor (As) as constraints 
(column 10). In this case the standard deviations for all the parameters are 
smaller (or at the most, equal) compared to those in any other column of Table 5, 

Further, it is also noticeable that the improvement from a non -constrained 
solution to a constrained solution, both with three or four constraints, is more 
significant compared to the improvement from a constrained solution using variances 
only to a constrained solution using (3 X 3) banded diagonal or full variance-covariance 
matrix. The improvement from the solution using (3 X 3) banded diagonal to the 
solution using full variance-covariance matrix is, however, marginal. 

A study of Table 6 indicates in all the three cases an overall improvement in 
correlation from a non-constrained to a constrained solution with four constraints 
(three rotations and one scale factor). The improvement in correlation between 
translations and rotations is quite significant while the same in not reflected between 
translations and scale factor. However, the improvement pattern from Table 7 is 
not straightforward. The correlations between translations and rotations show a 
downward trend from the solutions using variances only to the solutions using full 
variance -covariance matrix in all the three cases while the correlations between 
translations and As show an upward trend. 
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APPENDIX I 


Fortran IV Program with Subroutines 



♦♦♦♦♦♦♦♦ Ht 
♦♦♦♦♦♦♦♦ # 

39c3$c;4c3Cc3^i$c:4(4c ♦♦♦♦♦♦♦♦ # 

4c4c4c4c4c4c4i4t ♦♦♦♦♦♦♦♦ # 

4c4^4c^«4c4:« TRANSFORMATION OF AXES ♦♦♦♦♦♦♦♦ « 

4c***«*4r* ♦♦♦♦«♦♦ ^ 

4c4c« 4c 4c* *4c«<«<<'<'««4c4r4E:4^«3t^4c4c4c*«^4t4c3^*««*#«*«4c4c*^«^**4t4(#^««4^4c«««««4t« ♦♦♦♦♦♦ ♦ 

3^4c4(*4c4t4c*«««««*:fr***4c*4c********«*»*4(**4c4c«******A(*4t4c««**«««««**#%*#4c^ ♦ 


**4c***«* 
******** 
******** 
******** 
******** 
******** 
******** 
******** 
f ******** 

C »**^S<^**** 
******** 
******** 
******** 
******** 
******** 
******** 
******** 


PROGRAM WORKS ON TWO SETS OF COORDINATES — EITHER 
SET CAN BE INPUT AS ELLIPSOIDAL COORDINATES, BOTH IN 
DEGREES AND METERS OR IN GEOS FORMAT, IN SUCH CASE 
SEMI-MAJOR AXIS 'A* AND ECCENTRICITY •£• ARE NEEDED. 

UVW MATRIX TAKES COORDINATES IN THE FIRST SYSTEM 
f IN FORMAT I5,3F15.5 ) 

XYZ MATRIX TAKES COORDINATES IN THE SECOND SYSTEM 
( IN FORMAT I4,5X,OF16.5 ) 

MAXIMUM NUMBER OF INPUT POINTS FOR EACH SYSTEM 35 


******** * 
******** * 
»♦*♦ * 
**** * 
**** * 
**** * 
***♦ * 
*♦** * 
**** * 
**** * 
♦♦** ♦ 
**** * 
*♦♦♦ ♦ 
**** * 
**** * 
***♦ * 
******** * 


********«««««****♦*♦♦*♦*♦**♦♦♦♦*♦♦♦♦♦♦**♦***♦♦♦♦♦*♦*♦«««««♦♦♦♦♦♦★♦ ♦ 
********<<<<«««****************************************«««««*♦****** ♦ 
********<««<««*♦♦*****♦♦**♦*♦***♦****♦**♦**♦*♦********«««««♦♦♦♦*♦♦♦ ♦ 


SUBROUTINE 


•CSTRNT* 


******** * 
******** * 
**** * 
*♦♦* * 
*♦** * 
*♦** « 
♦*♦* * 
***♦ ♦ 
**** * 
**** * 
**** * 
***♦ * 
**** ♦ 
***♦ * 
**** * 
******** * 
******** * 

4c*** 4c 4c** «««««** *********************** ***************««<<<«<******** ♦ 

********«««««********♦***♦♦♦»♦**♦♦♦♦♦♦*♦»♦*♦♦***♦**♦*«««««♦,#***♦♦* ♦ 


******** 
******** 
******** 
{; ******** 
c *4****** 
******** 
******** 
******** 
******** 
******** 
4c 4t ****** 
******** 
******** 
*** ***** 
c ******** 
C ♦«*3^c**** 

C ******** 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


SOLVFS FOR TRANSFORMATION CASE WHEN CONSTRAINTS ARE 
TO BE APPLIED FOR THREE ROTATIONS. NECESSARY COUNTER 
~ KCOOE(ll) — IS TO BE COOED AS • 4 » . 


TWO SOLUTIONS ARE OBTAINED WITH THE SAME DATA 

FIRST WITHOUT CONSTRAINTS AND SECOND WITH CONSTRAINTS. 

INPUT CONSTRAINTS ARE OBTAINED FROM SUBROUTINE •EULERS* 
AND SUBROUTINE • SCALE*. 


SUBROUTINE 


•TFORM* 


******** * 
**** * 
**** ♦ 
**♦♦ » 
**** * 
**** * 
******** * 
******** * 

*4c ******««<««<♦«*****♦**♦*****♦*♦*****♦***♦*♦♦****♦♦**«««««♦♦***♦** * 

4c4c***4c4c*<<««««#^<c**»**<f**«******^^#**^<t<t^**«**»^+<t<t#**«««««******** * 


******** 
****4;* ** 
******** 
******** 
******** 
******** 
******** 
******** 


TRANSFORMATION PARAMETERS SOLVED UNDER THREE CASES. 
REFER KCOOEO) ALSO. 
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«4c<r4r4(#^« ♦ 

*♦*♦ « 
♦ ♦♦♦ ♦ 

♦ ♦«« ♦ 
❖ ♦♦♦ ♦ 

♦ ♦♦♦ ♦ 
♦♦♦♦ ♦ 


4c 4: 4^4: 3^ 4: 4: 4r 
4c 4^4^41^3^4:4 
4r4f4«4444 
444*44*4 
44444444 
44444444 
44444444 
44444444 
44444444 
44444444 
44444444 
44444444 
44444444 
44444444 


SUBROUTIve — 'SCALF* 


SCALF FACTOR BFTWFFN SYSTEM #1 AMD SYSTEM #2 IS 
COMPUTED BY COMPARISON OF U-K)TH CHORDS IN THE 
TWO SYSTEMS. THE WEIGHT FOR EACH ESTIMATE OF 
SCALE FACTOR CCRRESPOND ING TO (I-K)TH CHORD IS 
COMPUTED USING VARIANCE - COVARIANCE MATRICES 
OF •!• TH AND »K« TH POINTS OF SYSTEM #1 AND 
SYSTEM ^2 THROUGH ERROR PROROGATION. 


44 444444 <<<<<<<<<<44444444 4444444444 4444444444444444444444««<««<4«44*44* ♦ 
44444444<<<<<<<<<<4444444444444444444444444444444444444444<<<<<<<<<<44444444 4 


444:«34444 
444A4444 
44444444 
44444444 
4 4444444 
44444444 
(; 44444444 

44444444 

44444444 
44444444 
44444444 
44444444 


SUBROUTINE ~ •EULERS* 


EULERS ANGLES ARE COMPUTED FROM DIRECTION COSINES. 
WEIGHT MATRIX 'P* FOR »TIK AND OIK« IS COMPUTED 
USING VARIANCE - COVARIANCE MATRICES OF THE POINT 
COORDINATES OF THE SYSTEMS BY ERROR PROROGATION . 


*4*4444* 4 

44444444 * 

44444444 * 
444* ♦ 

444« « 
44*4 ♦ 
«44« « 

4444 4 
4444 4 
44444444 * 

44444444 * 


4 4n3 44 444«««<«<44 4 4 4*44 44 44444444444*444^444444444*4444<««««<4*«***## ♦ 
44444 444 <<«««< <4444 4444 *444444444*444444*4 44444444*4444«««««44*«^^## ♦ 


44444444 
44444444 
44444444 
44444444 

44444444 

44444444 
44444444 
44444444 
44444444 
44444444 
44444444 


SUBROUTINE 


•SETUP' 


SETS UP AND STORES WEIGHT MATRIX 'P' IN VECTOR FORM. 


44444444 * 

44444444 * 

4444 ♦ 
4444 4 
4444 4 
4444 4 
44444444 4 
44444444 4 
44444444 4 
44444444 4 
44444444 4 

4444444V<<<<<<<<<<4444444444444444444444444444444444444444<<<<<<<<<<44444444 4 
44444444<<<<<<<<<<4444444444444444444444444444444444444444<<<<<<<<<<44444444 4 
C 44*44444 44444444 4 

C 44444444 44444444 4 

44444444 *44*4444 * 

44444444 INPUT KCOOES 44444444 4 

44444444 JO BF CODED WITH EACH DATASET 44444444 4 

44444444 44444444 4 

44444444 44444*44 * 

C 44444444 44444444 4 

C 44444444 4*444444 4 

C 44444444 444*4444 4 

C 4444 4444«««<«<*******4**********4**4*****4******44**4*«««««44*444** 4 

C 44 444 444<<««««******4* ***************** *44****4**4****^<^^^<^<^^^^^1^^^^*^^ * 
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C <<<<<<<<<< ♦♦«♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦»♦««««<<♦♦♦♦★♦♦♦ ♦ 


c 





44444444 

4 

c 

5#cijr#:er3{i:> 




44444444 

4 

c 

:<r# 




44444444 

4 

c 





44444444 

4 

c 


KCODF^( 

1 ) 

= 'TaTAL NUMBER OF POINTS — IN (12) FIELD. 

4444 

4 

c 





4444 

4 

c 


KCODE ( 

2 ) 

= 'PAPAMETERS RECUIPED IN THE SOLUTION* 

4444 

4 

c 




3 DENOTES ONLY TRANSLATIONS OR 

4444 

4 

c 




ROTATIONS — SEE KC00E(14) 

4444 

4 

c 




TO SEPARATE THESE SOLUTIONS. 

4444 

4 

c 

4c 3^ 4 c34c 



4 DENOTES TRANSLATIONS AND SCALE. 

4444 

4 

c 

4t# 



7 DENOTES TRANSLATION, SCALE AND 

4444 

4 

c 

:4^4M>s^*4:3{r4t 



THREE ROTATIONS. 

4444 

4 

c 

4c*4:*3Ct4^4c3(c 




4444 

4 

c 

:^4c 4(4r^4r3^4[ 

KCOOE ( 

3) 

= 'WHETHER CONSTRAINED SOLUTION IS NFFDED* 

4444 

4 

c 

4f4t4t4r4fJjs4:5X 



0 DENOTES NO SUCH SOLUTION 

4444 

4 

c 

*4c4i 4f:«f «:3fC3^ 



1 OFNOTFS CONSTRAINFD SOLUTION. 

4444 

4 

c 

4c4i;^:3Cc4t4:** 




4444 

4 

c 

4i 4c4: 

KCODE ( 

4) 

= 'FIRST SYSTEM IN ELLIPSOIDAL COORDINATES 

4444 

4 

c 

*♦ :^<c;^4:4t4t 



IN DEGREES FOP PHI AND LEMDA — 

4444 

4 

c 




AND METERS FOR HEIGHTS* 

4444 

4 

c 

4c4c4c4(’^^^>}( 



0 DENOTES NO SUCH CASE 

4444 

4 

c 

4c 4: 4. 4i 4(3^ 4c 4c 



1 DENOTES SUCH INPUT 

4444 

4 

c 

4c4c4.4c4c4^3«^2#r 




4444 

4 

c 

4c 4c 4; sjf 4 :j> 

KCODE ( 

5) 

= 'INPUT FOR FIRST SYSTEM IN GEOS FORMAT 

4444 

4 

c 

«4c 4c 44c 4c 3^4 



0 DENOTES NO SUCH CASE 

4444 

4 

c 

444:^:50t:t-^3fr 



! DENOTES SUCH INPUT 

4444 

4 

c 

44 444t 4 44 




4444 

4 

c 

444444 

KCODE ( 

6 ) 

= 'SECOND SYSTEM IN ELLIPSOIDAL COORDINATES 

4444 

4 

c 

44444444 



IN DEGREES FOR PHI AND LEMDA — 

4444 

4 

c 

444444 44 



AND HEIGHTS IN METERS* 

4444 

4 

c 

44444444 



0 DENOTES NO SUCH CASE 

4444 

4 

c 

44444444 



1 DENOTES SUCH INPUT 

4444 

4 

c 

44444444 




4444 

4 

c 

44444444 

KCODE ( 

7) 

,ir.. • input FOR SECOND SYSTEM IN (iEOS FORMAT* 

4444 

4 

c 

44444444 



0 DENOTES NO SUCH CASE 

4444 

4 

c 

44444444 



1 DENOTES SUCH INPUT 

4444 

4 

c 

44444444 




4444 

4 

c 

44444444 

KCODE ( 

8 ) 

= 'VARIANCE COVARIANCE MATRIX AS DIAGONAL* 

4444 

4 

c 

44444444 



IN FORMAT *3F10.2 — 1 CARDS PER STN. 

4444 

4 

r 

44444444 



0 DENOTES NO SUCH CASE 

4444 

4 

c 

44444444 



1 DENOTES SUCH INPUT 

4444 

4 

c 

44444444 




4444 

4 

c 

44444444 

KCDDE( 

9) 

= 'VARIANCE - COVARIANCE MATRIX IN 3X3 FORM* 

4444 

4 

c 

44444444 



0 DENOTES NO SUCH CASE 

4444 

4 

c 

444 44444 



IN FORMAT *3F5.2* — 3 CARDS PER STN. 

4444 

4 

c 

44444444 



1 DENOTES SUCH INPUT 

4444 

4 

c 

44444444 




4444 

4 

c 

44444444 

KCODE ( 10 ) 

= 'VARIANCE - COVARIANCE MATRIX IN FULL AS 

4444 

4 

c 

44444444 



UPPER TRIANGLE (ROW-WISE) IN VECTOR FORM* 

4444 

4 

c 

44444444 



IN FORMAT *8FI0.4* — EACH NEW ROW TO 

4444 

4 

c 

44444444 



BEGIN ON A NEW CARD FROM COLUMN 1. 

4444 

4 

c 

44444444 



SEE KCODE( 12) ALSO. 

4444 

4 

c 

44444444 



0 DENOTES NO SUCH CASE 

4444 

4 

c 

44444444 



1 DENOTES SUCH INPUTSE 

4444 

4 

c 

444444*4 




4444 

4 
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«:^«: Qc 4 cA 3$^4 

!0c<c4^^4r4c4c^ 

4r4c«4c:^:^a9r« 

4c4t44t««<c* 

4:^«]ft4c«3^4c 

4r4c4 4:^4c^:4r 
Xc:^i4c30:i^i(r«^ 

«3^t««««l}t^t 


KCODE(ll) = ®FO« CONSTRAINED SOLUTION TO BE CODED A$o 
»3® ROTATIONS ARE CONSTRAINED o 

ROTATIONS AND SCALE ARE CONSTRAINED 

KCODE(12) = ®TO BE CODED »1» TO OMIT CORRELATION 

WITH KC0DEI101 AS AN ALTERNATE SOLUTIONo 


KCODE(13) = ® 


TOTAL NUMBER OF TRANSFORMATIONS 

TO BE PERFORMED ® TO BE CODED 

WITH THE LAST DATA SET IN (12) FIELOo 


KCODE(l^) 


®3 PARAMETER SOLUTION ONLY® 


SOLUTION FOR TRANSLATIONS 
SOLUTION FOR ROTATIONS 


0000 o 
0000 o 
^ oooo o 
0000 o 
0000 o 
oooo ❖ 
oooo ♦ 
oooo o 
oooo o 
oooo o 
oooo o 
oooo o 
<I000 o 
oooo 
oooo 
oooo 
oooooooo 
oooooooo 
oooooooo 


<SOOOOOOOOOOOOOOOOOOOOOOOO<<<<<<<<<<OOOOOOOO 
00000000000000000000000000000000<<<<<<<<<<00000000 
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IMPLICIT REAL * 8(A-H , 0-Z) 

PEAL *8 LEMDA,NI,M02 

n I MEN SI ON XYZ(35,3) fRANGLE(^) ,VR0T(A,4) ,NAHEl(3) , 

2 A (3600) ,W( 1200) ,P( 2400) ,UVW( 35,3) ,NAME2 (3 ) , 

3 AA(3,105) ,BEi(3,105) ,NSTA(35) ,KSTA{35) ,KC0DE( 15) 

COMMON /WEIGHT/ P 

COMMON /CODE/ KCOOE 

COMMON /INAME/ NAME1,NAHE2 

COMMON NSTA,KSTA,NN,NM,UVW,XYZ , A,W,KPR,KPARM 

COMMON /ANGLE/ RANGLE,VROT 

DATA MINUS/IH-/ 

PIl = 3.141592653589793DC 

PHO = 180. DO/PI I 

RHOS = RHO+3600.D0 

KDUNT = 1 


♦ ♦****=jt* t?FAn IN VARIOUS CODES INVOLVED 


COO READ (5, 1) (KCODE(I), I = 1 , 15 ) , ( NAMEl ( 1 ) , I =1 , 3 ) , 

2 (NAMF2(I),I=1,3) 

1 FORMAT ( I2,11I1,I2,2I1,3X,3A4,3X,3A4) 

WRITE (6, 2) (KCQDE(I), I = 1,15) 

2 FORMAT ( 'I', //////, 25X, •KCOOE INPUT' , //, 20X , 15 12 , // ) 

NO = KCODE(l) 

IP (KC0DE(4) .EO.O.AND.KCODE(5 ) .EQ.O) GO TO 12 

<.****««;» read IN DATA FOR THE FIRST SYSTEM 

RPAO (5, 3) AE1,F 

3 FORMAT (2F15.10) 

F = l.DO/F 

E2 = 2.D04F - F*F 

IF ( KCnnE(5) ,E0. 1) GO TO 6 

read in ELLIPSOIDAL COORDINATES IN DEGREES AND HEIGHT 
DO 5 I = 1 , NO 

READ (5, 4) NSTA( I) ,PHI,LEMDA,HT 

4 FORMAT ( I4,5X,3F16.9) 

PHI = PHI / RHO 

LEMDA = LEMDA / RHO 

WW = (1.D0-E2 *DSIN( PHI )«DSIN(PHI ) )**0.5DO 

UVW(I,1)= (AF1/WW+HT)>ADC0S(PHI)*DC0S(LEM0A) 

UVW( I ,2) = (AE1/WW^HT)>X0C0S(PHI )*DS ini LEMDA) 

UVW(I,3)= ( ( { AE1*(1.D0-E2 ) ) /WW) +HT ) A-DSIN ( PHI ) 

5 continue 

GO TO 15 

C 

C 

E ♦*♦#«**« read in ellipsoidal COORDINATES IN GEOS FORMAT 

C 

C 

6 DO 11 I = 1 , NO 

READ (5 , 7) NSTAd ),ISN,IPH,MPH,SPH,ILM,MLM,SLM,HT 

7 FORMAT ( I4,20X,A1,2I3,F8.3,213,F8.3,F10.2) 

LEMDA = ( 1LM+ ( (MLM+( SLM/60.D0) )/60.D0) )/RHO 


- 32 - 



on ooono ooooo oonoo 


IF (ISN .EC. MINUS) GO TO 8 

PHI = ( IPH+( (MPH*(SPH/60.00) )/60.00) )/PHO 

GO TO 10 

8 PHI = -IIPH+IIHPH^JSPH/60.DO))/60.00))/RHO 

10 WW = (l.DO-EZ^nSINIPHI )*nSlN(PHI) )**0.5D0 

UVW(Itl) = (AEl/WW+HT)*nCOS<PHI)POCOSaEMOA) 

UVW(I.2) = (AEl/WW*HT)*nCOS(PHI)*DSIN(LEMOA) 

UVW(I*3) = ( ( (AE1*(1.D0-E2) )/WW)-»HT)*DSIN(PHI) 

11 CONTINUE 
GO TO 15 

C 

C 

c *♦**«♦** read in rectangular coordinates ( Uf Vt w ) in meters 

c 

c 

12 DO 14 I = 1 , NO 

READ (5, 13) NSTA(I),lUVW)I,J),Jsl,3) 

13 FORMAT! I4,5Xf3F16.5) 

14 CONTINUE 


READ IN COORDINATES OF THE SECOND SYSTEM 


15 IF (KC0DE(6) .EO.l.OR .KCOOE (7) .EO. 1 ) GO TO 20 


*««««««« read in RECTANGULAR COORDINATES ( X« Yt Z ) IN METERS 


16 

DO 18 I = 

RFAD 

FORMAT 

1 f NO 

C5, 16) KSTAm,(XY2(I,J), J=l,3) 
(lA»5Xf3F16.9) 

18 

CONTINUE 


20 

GO TO AO 
READ 

(5» 22) AE2tF 

22 

FORMAT 

(2F15*10) 


F = 

E2 

IF 

l.DO/F 

2,D0«F - F*F 

( KC0DE(7) .EQ. 1) GO TO 25 


«««««**« r^AD in ELLIPSOIDAL COORDINATES IN DEGREES AND HEIGHT 


DO 24 I = 1 , NO 

READ (5, 23) KSTAI 1 ) ,PH 1 , LFMOA ,HT 

23 FORMAT ( 14 , 5X t 3F16 .9 ) 

PHI = PHI / RHO 

LEMDA = LEMDA / RHO 

HW = (1.D0-E2 •DSlN(PHI)*nSIN(PHI))»*0.5D0 

XYZIIf 1)= ( AE2/WW+HT » *DCOS( PHI) *OCOS( LEMDA) 

XYZIIt2)= (AE2/WH+HT)’»0CCSIPH1 )*DS1N(LEMDA) . 
XYZ(I,3)= mAE2-M1.00-E2 ) ) /HW) fHTIPDSlNI PHI ) 

24 CONTINUE 

GO TO 40 
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c RCio IN FLL1PSDID4L CDOROINATES IN GEOS FORMAT 

C 

C 

25 no 31 1 = 1 * NO 

READ (5, 26) KST A ( I ) , I SN, IPH.MPHf SPH , I LM »MLM t SLM |HT 

26 FORMAT ( 1 4 , 20X , A1 , 213 , F8 .3 f 2 1 3 , F8 . 3 t F 10.2 ) 

LEMCA = ( ILM-f { (MLM+ISLM/60.00) )/60.D0) )/RHO 

IF (ISN .EO. MINUS) GO TO 28 

PHI = (IPH+I (MPH+(SPH/60.D0) )/60.D0) )/RHO 

GO TO 30 

28 PHI = -( IPH'fl (MPH+(SPH/60.D0) )/60.D0) )/RHO 

30 WW = (1.00-E2*nSIN(PHI)*OSlNfPHI))**O.5D0 

XY7(1,1) = ( AF2/WW+HT)*nC0S(PHI)*DC0S(LEM0A) 

XYZ(1,2) = (AF2/WW+HT)*DC0S(PHI)*DSIN(LEMDA) 

XYZtI,3) = ( ( (AE2*ll.D0-F2))/WW)+HT)*DSIN(PHI) 

31 CONTINUE 
C 

C 

C 

C V^PITING OF RFAD IN DATA FOR THF TWO SYSTEM IN RECTANGULAR COORDINATES 

C 

C 

r 

40 WP, ITF(6, 42) 

42 FORMAT! 'IS///. 25X, 'RECTANGULAR COORDINATES FOR FIRST SYSTEM',///) 
WPITr(6, 43) 

43 FORMAT!' ' , 1 3 X , ' STN.NO . ' , 1 2X , ' U ' , 13X , ' V ' , 16X , ' W' ,/ ) 

no 46 I = 1 , NO 

W5ITF16, 44) NSTA!1), !UVW!I,J), J=l,3) 

44 FORMAT!' • ,13X,I5,F20.4,2F16.4, !14X,I5,F20.4,2F16.4)) 

46 CONtTINUF 

W'>1TR!6,50) 

50 FORMAT! ' I' ,///,25X, 'RECTANGULAR COORDINATES FOR SECOND SYSTEM',/) 
WRIT£!6,52) 

52 FORMAT!'-' , 1 3X ,' STN.NO 12X ,' X ', 13X , 'Y ', 16X ,' Z ',/ ) 

DO 60 1 = 1 , NO 

WRITF!6, 58) KSTAIU, !XYZ|1,J), J=l,3) 

58 FORMAT!' ',13X,I5,F20.4,2F16.4,!14X,I5,F2D.4,2F16.4)) 

60 CONTINUF 
C 
r 

C SEPARATING THE TYPE OF SOLUTION REQUIRED 

C 

c 

KPARM = KCODE(ll) 


IF 

(KC0DF(8) 

.NF. 1) 

GO 

TO 

62 

KPR 

= 1 





GO TO 

75 





62 IE 

(KC0DE(9) 

.NF. 1) 

GO 

TO 

64 


= 2 

GC TO 75 
64 KPR = 3 

IF (KCDDE(10).EQ.1.AND.KCODEI12).E0.1) KPR = 2 

75 MM = NO “ 1 

m = NO ♦ NM 

NNM = 3«N0 

IF (KC0DE(14) .EQ. 0) GO TO 85 

CALL PULERS ( NO ,NNN , AA , BB ) 
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IF (KC0DE(14) .E0.1.AND.KC0DEI2) .E0.3) GO TO 95 

85 CALL TFORM (NO,NNN) 

IF (KOUNT .EQ. KCODEI13)) GO TO 95 
KOUNT = KOUNT ♦ 1 
GO TO 1000 
95 STOP 
ENO 
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C 

C 


C ♦ 

c 'EULFR'S ANGLES* ♦♦♦ ♦ 

C * 

C * 

f :(t#*«*:}f+««*#«**»««**«^**^*###:>j4t*^^*****#*^^*«:0t*****^*^*##*^******^***^^«*»**** 

C ♦ 

c ♦ 

c RETWEEN TWO COORDINATES SYSTEMS COMPUTED FROM OIRFCTTON COSINES ♦♦♦♦ ♦ 

C ♦ 

C ♦ 

C ♦ 


C 

C 

C 

C 

C 

C 

SUBROUTINE EULERS ( NO , NNN» AA » BB ) 

IMPLICIT REAL ♦S CA-H,0-Z) 

REAL *8 NItN7,M02 

DIMENSION UVW(35,3),XYZ(35»3 ) * A ( 3600 ) ,W ( 1200 ), NAME 1 C3 ) f 
2P1(6,6) »G(2,6) ,GP(2,6) tGT(6,2) tPP(2»2) fKX(2) ,KYC2) ,NAME2(3) f 
3B(2,4) ,BT( 4,2) ,P2( 6,6) , INDEX (40) , I NV (40 ) ,OXYZ( 4500 ) ,NZ ( 4,4) , 
4P(2400),BS(2,4),KSTA(35) , NSTA ( 35 ) , 0UVW( 4500 ) ,AA(3,NNN) ,BB(3,NNN), 
5P0(2,2) ,PR(4,4),NI(3,3),nX(3),U(3),VAR(3,3),K0(3),KC00E(15),L0C3) 
COMMON /WEIGHT/ P 

COMMON /CODE/ KCCDF 

COMMON /ANGLE/ S,DX,NZ 

COMMON /INAME/ NAME1,NAME2 

COMMON /SEAC/ DW ,D S , DA 1 , DB 1 , DC 1 , 

2 nA2,DB2,DC2,RIKl ,RIK2,P1,P2 

COMMON NSTA,KSTA,NN,NM,UVW,XYZ,A,W,KPR,KPARM 

PI I = 3»14159265358P793D0 

RHO = 18C.D0/PII 

RHOS = RHO*3600.DO 

nw = 0.00 

DS = O.DO 

S = O.DO 

VSF = 0.00 

WT = 0.00 

LL = 1 


d S’^TTING UP Of MATRIX »B* — COMMON TO ALL SOLUTION 

C 

C 

c 


B { 1 , 1 ) 

= “l.DO 

B ( 1 . 2 ) 

= O.DO 

8(1.3) 

= 1.00 

P(l,^) 

= O.DO 

R(2,l) 

= O.DO 

R(2,2) 

II 

1 

• 

a 

o 
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= O.DO 

NJS = m/2 

B(2,4) = 1.00 

DC 1 I = 1 f 2 

nn 1 J = 1 T 4 

PT(Jtl) =r B(I,J) 

1 CCNTINUF 

DO 2 1 = 1 ♦ A 

DC 2 J = 1,4 

PR(1,J) = O.DO 
2 COMTINUF 

IF (KC0DE(8) .E0.1.0R,KCnDE(9) .EC.l ) GO TO 10 

C 
C 
C 


C ♦ 

c ♦ 

C FULL VAR IANCE^COVARIANCE CASE ♦♦♦♦ ♦ 

C ♦ 

c ♦ 


C 

C 

c 

C ^FADING IN VARIANCE--CDVARIANCES FOR 'FIRST SYSTEM* 

C 

C 

C 

JK = 1 

f?C) 6 I = 1 , NNN 

JL = JK ♦ NNN - I 

RFAD (5, 3) (OUVW(J), J = JK,JL) 

3 FORMAT (8F10.4) 

or 4 L = LL » 3 
P1(LL,L)= QUVW( JK+L-LL) 

4 P1(L,LL)= P1(LL,L) 

WRITE m (PKLLtM) ♦ M = 1 , 3) 

LL = LL ♦ 1 

IF (LL .EQ. 4) LL = 1 

6 JK = JL ♦ 1 

PFWINO 1 

C 

C 

C pfAOING IN VARIANCE-COVARIANCES FOR *SECOND SYSTEM* 

C 

C 

c 

LL = 1 

JK = 1 

DO 9 I = 1 , NNN 

JL = JK ♦ NNN - I 

READ (5, 7) (OXYZ(J), J = JK,JL) 

7 FORMAT (8F10.4) 

DC B L = LL , 3 
P?(LL,L)= OXYZ( JK+L-LL) 

8 R?(L»LL)= P2(LL,L) 

WRITE (2) (P2(LL»M), M = 1 , 3) 
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LL = LL ♦ 1 

(LL .FO. 4) LL = 1 
9 JK = JL ♦ 1 

^^FWIND 2 

GC TC 24 
C 
C 


c ♦ 

c ♦ 

C DIAGONAL 0=^ 3X3 BANDED CASE ♦♦♦♦ ♦ 

C ♦ 

C * 


C 

C 

C 

Q rfading in variancf-covariancf for first system 

c 

c 

10 Dn 17 I = 1 t NO 

KK = ( I-n*3 + 1 

KM = KK + 2 

IF (KCnOF(e) .FO. 1) GO TO 13 

c 

c 

c 

c VARIANCE - COVARIANCE MATRIX IN 3X3 BANDED FORM 

C 

c 

c 

c 

00 12 J = 1*3 

READ (5,11) (AA(J,K), K = KK*KM) 

11 Ff?RMAT(3F5.2) 

12 WRITF(l) (AA(J,K), K=:KK,KM) 

GO TO 17 

C 

c 

(;^ VARIANCE - COVARIANCE MATRIX IN DIAGONAL FORM (ONLY VARIANCES) 

C 

C 

13 00 14 J = 1 , 3 
no 14 K = KK * KM 

14 AA(J,K) = C.DO 

read (5,15) (AA(K,(K + KK--1)), K = 1,3) 

15 FORMAT (3F10.2) 

DO 16 J = 1 , 3 

16 Wf^lTF(l) (AA(J,K), K=KK,KM) 

17 CONTINUE 
REWIND I 

C 

C 

Q READING IN VARl ANCF-COVARI ANCE FOR SECOND SYSTEM 

C 

C 

DC ?3 1 = 1 , NO 

KK = ( i-j)^3 + 1 
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KM = KK 2 

IP (KC0DF(8) .EO. 1) GO TO 20 


♦ *'** VARIANCE - COVARIANCE MATRIX IN 3X3 BANDED FORM 


Df,! 19 J = 1 , 3 

RFAD (5fl8) K=KK,KM> 

16 FORMAT (3F5.2) 

19 WP1TE(2) (PRU*K), K=KK,KM) 

GO TO 23 


VARIANCE - COVARIANCE MATRIX IN DIAGONAL FORM (ONLY VARIANCES) 


20 DO 21 J = 1 , 3 
DC 21 K = KK,KM 

21 PB(J,K) = O.DO 

READ (5,15) (BB(K,(K+KK-l) ), K = 1,3) 

on 22 J = 1 , 3 

22 WPITF(2) (BB(J,K), K=KK,KM) 

23 continue 
REWIND 2 

C 

c 


C * 

C >»*** ECRMING MATRICES 'A*, »W', AND 'P' FOR THE ENTIRE SYSTEM * 

c by computing direction cosines for each line betneen *♦** * 

c ***♦ ANY ONE SET OF TWO GIVEN POINTS. **** * 

C * 

C • 


c 

c 

74 MKR = 1 

KMT = 1 

MK = I 

INDfcXU) = I 

MMl = NNN + 1 

on 75 I = 1 , NP 

25’ INV( I ) = 3«I - 1 

nr 5C I = 1 , NM 

DC 26 J = 1 , 6 

on 26 K = 1 # 6 

P1(J»K) = O.DC 
P2(J,K) = O.DO 
26 CrNTINUP 

IF (KCODEdO.FO. 1) GC TO 28 

nn 27 J = 1*3 

DO 27 L = 1*3 

LL = (1-1) # 3 + L 

PlUtL) = AA(J*LL) 

P2(J*L) = BB(J*LU 
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?7 CCMTINUE 

GO TO 32 

ZB LL = INDEX(I) 

no 30 J = 1 f 3 

DO 29 L = J t 3 

LLL = LL + L - J 

PI { J»L) = CUVW(LLL) 

2^ P2(JtL) = QXYZ(LLL) 

MMl = MMl - 1 

30 LL = LL + MMl 

32 JJ =1 + 1 

INOFX(JJ)= LL 
MM? = MMl 

no 50 K = JJ ♦ NO 

IF <KC0DE(8 >.FQ,1 .OH . KCODF f 9 ) . EQ . 1 > GO TO ^3 

LL = INDEX(K) 

on 3A J = 4 ♦ 6 

DO 33 L = J , 6 

LLL = LL + L - J 

Pl(JfL) = QUVW(LLL) 

53 P?(J,L) = OXYZ(LLL) 

MM? = MM2 - 1 

34 LL = LL + MM2 

KP = K + 1 

1ND5X(KP)= LL 

III = INOEXd) ♦ INV(K^I) 

IF (KCC0E(12) .FQ. 1) GO TO 41 

nn 3P J = 1 » 3 

DO 36 L = 4 t 6 

LLL = III + L - 3 

PKJtL) = QUVW(LLL) 

36 P?(J,L) = OXYZILLL) 

38 III = III + (NNN -(3*( I-l) )-J) 

41 nn 42 J = 1,6 

nr 42 L = 1,6 

Pl(LtJ) = P1(J,L) 

4? P?(Lf J) - P2(J,L) 

GO TO 45 

43 DO 44 L = 4,6 

JKL = L 3 

or 44 M = 4,6 

KLM = {K-2)^3 + M 

P1(L,M) = AA(JKL,KLM) 

= BBfJKL,KLM) 

44 CONTINUE 

45 KFM = MKR + NN 

KMS = MKR + (2*NN) 

C 

0 

c COMPUTING DIRECTION COSINES FOR FIRST SYSTEM 

C 

C 

DAI = UVW(K,1 ) - UVW( 1,1) 

PB1 = UVW(K,2) -- UVW(I,2) 

DCl - UVW(K,3) - UVW(I,3) 

oiKl = n$QRT(DAl*DAl+DBl«DBl+DCl^DCl) 

AIKI = DAI/RIKI 
PlKl = DBl/RIKl 
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CIKI = DCl/RIKl 

TIKI = -0ATAN2(B1K1,AIK1) 

IF (TIKl.LT.O.) TIKI =C360,00+TIK1^RHO) /RHO 

ABl = DS0RT(AIK1^AIK1+BIK1»BIK1) 

DIKI = OATAN2(CIKltABn 


COMPUTING ni'RFGTION COSINFS FOR SFCOND SYSTFM 


DA 2 

OB 2 

DC2 

RIK2 

AIK2 

B1K2 

CIK2 

TIK? 

IF 

AP? 

DIK.2 


♦♦♦♦ PETTING UP MATRICES 'A* AND ’H* — COMMON TO ALL SOLUTION 


A(MKR) = l.DO 
A(MKR+1)= O.DO 

A(KSM) = DSIN(TIK2I*0TAN(0IK2> 
A(KSM-m= 0C0S(TIK2) 

A(KMS) = -DC0S(TIK2)^0TAN(DIK2) 
A(KMS + 1)= nSIMTIK2) 

W(MKR) = TIKI - TIK2 
DIKI “ DIK2 
C 


C ♦ 
C ♦ 

C ♦♦♦♦ forming VAR-COVARIANCE MATRIX FOR ‘TIK* AND •DIK» ♦♦♦♦ ♦ 
C ♦♦♦« THROUGH PROROGATION OF ERRORS — WHERE »TIK» AND ♦♦♦♦ ♦ 
C ARE GEODETIC HOUR ANGLE AND DECLINATION. ♦♦♦♦ ♦ 
C ♦ 
C ♦ 


C 

C 

C 

C FIRST SYSTEM ♦♦♦♦ ♦ 

C 

C 

C 

C 

OABl = DA1«0A140B1*DB1 

DBA = OSQRT(OABl) 

G<1,11 = -OBl/DABl 

G(l,2) = DAl/DABl 

G(l,3) = O.DO 

r,n,4) = -Ga,i) 
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= XYZ(Kfl) - XYZn*l) 

= XYZ(K,2) - XYZ(If2) 

^ XYZ(K,3) - XYZ(If3l 
= DSORT(OA2*DA2^DB2^DB2+DC2«OC2) 

= DA2/RIK2 
= DB2/RIK2 
= DC2/R1K2 
= -DATAN2(BIK?,AIK2) 

(TIK2.lt. 0. ) TIK2 = (360 .DO^TI K 2^RH0) /RHO 
= DSQRT( AIK2»AIK2+BIK2*BIK2) 

= DATAN2(ClK2tAB2l 



n o o o o o 


C 

C 

C 

C 

c 

c 

c 


^6 




A7 


G(l,5) 

-G(lt2) 

G n , 6 ) = 

O.DO 

GiZfl) = 

DA1*DC1/(DBA^P.IKI*PIK1) 

Gi?f 2) = 

Dei^OCl/(DBA«RIKl«RIKl) 

G ( 2 » 3 ) = 

-DBA/(RIK1*R1K1) 

G ( 2 f A ) = 

-G(2,l) 

G(2f5) = 

-G(2,2) 

G(?»6) = 

-G(2,3) 

DO A6 L = 

1 * 2 

DQ A6 M = 

1 » 6 

GT(M»L) = 

G(LrM) 

CONTINUE 


CALL 

DGMPRO(GTPl»GPt2t6,6) 

CALL 

OGMPRDCGP,GT,PP,2,6,2) 


SECOND 


DAP2 = 

0A2*DA2-*'0B2*DB2 

DAP = 

DSORT(DAB2) 

G(l,l) = 

-DB2/DAB2 

= 

nA2/DAB2 

G ( 1 , 3 ) = 

0.00 

G (1 , ) = 

-G(l,l) 

G(l,5) = 

~G(1,2) 

G ( 1 f IS ) = 

O.DO 

G(2,l) = 

DA2+DC2 / (DAB*RIK2*RIK2) 

G { 2 , 2 ) = 

DB2^DC2 / (DAR*RIK2*RIK2) 

G(2,3) = 

-DAB/CRIK2+RIK2) 

G(2.,4) = 

-G(2,l) 

G(2t5) = 

*G<2,2) 

G(2.ft) = 

«G(2,3) 

DO 47 L = 

1 f 2 

DO 47 M = 

1 » 6 

GT(M,L) = 

G(L>M) 

CCNTINUP 


CAUL 

DGMPRD(G»P2»GP»2,6,6) 

CALL 

DGMPRD(GP,GT,PQ»2f 6,2) 


SYSTEM 


♦ 


FORMING MATRIX »MI» FOR THE COMBINED SYSTEM 


no A8 L = 

1 f 2 

J 

L + 2 

nn AP M = 

1 • 2 

N 

M + 2 

PP. (LfM) = 

PQ(L,M) 

PR( JfN) = 

PP(L,M) 

CONTINUE 


CALL 

DGMPRD(B,PR,BS,2,A,A) 

CALL 

DGMPRDI BS,BT,PP,2,A,2) 

CALL 

DMINV(PP,2,DT,KX,KY) 
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P(KMT) = PP(1,1) 
P<KMT+} )= PP(2f n 
o(KMT^2 )= DP( lt2) 
P(KMT+3)= PPiZf?) 
MKR = MKR ♦ 2 
KMT = KHT ♦ 4 


IF (KCCOFUl) .FQ, 3) GO TO 50 

CALL SCALE ( IMS ,MK , S ,VSF ,WT) 

MK = MK ♦ 1 

C FINDING WFIGHTFD MEAN AND VARIANCE FOR ♦♦♦♦♦♦♦♦ « 

C *SCALF FACTOR* 6Y COMPARISON OF CHORDS IN ♦♦♦♦♦♦♦♦ « 

C THE TWO SYSTEMS BY CALLING SUBROUTINE •SCALE*. ♦♦♦♦♦♦♦♦ « 

C «'«««♦♦♦♦♦*♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦»♦♦ ♦♦♦♦♦♦<<««««♦♦♦♦♦♦♦♦ ♦ 

C 4r 4r 4c :(c ^ ♦ 

50 CONTINUE 


VSF 

0(i 75 J = 

IJK 

JKL = 

75 WRITPO) 
^tWlNO 
WRITFtA) 
PFWINO 


VSF ♦ 10. Dll 
1 » 3 

(J-1)*NN ♦ 1 
JJK ♦ NN - 1 
( A( I) , I s IJK 
3 

(W(K) , K=1,NN) 


JKL) 


FOPMING MATRIX 'N» AND INVFRTING THE SAME 



no Ro I = 

1 » 3 


RFAD (3) 

<W( J) ,J=ltNN) 


K1 

( I-1)«NN *► 1 


K2 . 

K1 + NN ^ 1 


MMM = 

0 


no 76 K := 

K1 , K2 


A(K) = 

O.DO 


LI 

( ((K-Kl)/2)+2) > 1 


L2 

LI + 1 


nn 7B L = 

LI , 12 


MMM = 

MMM + 1 

78 

A(K) 

A(K) + W(L)*PIMMM) 

FO 

CONTINUE 



DO 84 1 = 

1 f 4 


DO 84 J = 

1 t 4 


NZUtJ) = 

O.DO 


REWIND 

3 


DO 88 I = 

1 t 3 


P F AD ( 3 ) 

(W(L)t L=1 ♦ NN) 
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Pf/ 8f. J = 

1 ♦ 3 


N n J f 1) = 

O.DO 


DO 8*^ K = 

1 , NN 


III 

(J-D4NN ♦ K 

or 

N1 (Jt I ) = 

NI(J,n ♦ A(IU)^W(K) 

68 

CONTINUE 



REWIND 

3 


IF 

(KCODEdl ) .EO. 3) 


N7(l,l) = 

WT 

S9 

DP 91 1 = 

2 t A 


DO 91 J 

= ? f ^ 

91 

NZ(I,J) = 

NI(I-1,J-I) 


CALL 0!^INV(MI ,3,DF7»KQ,10) 


^ 4 ^ 4^ 4c 4c 3ft « «« 4c « « 3^ « 4 30c ♦ « 4c 4c ♦ * ♦ « « ♦ ♦ « « «« « 

♦ 

4c#*4c COMPUTING SOLUTION VECTOR • DX • FOR 3 ROTATION PARAMETERS ♦♦♦♦ ♦ 

A 

4 4c 4c 4^ 4c 3^ 4: 3^ 4r 4c 4c 3^ 3^ 3»: 4c4t « « 4( 3» 3tc 4c « 3ft :ft 4c 4r 3ft 4 « 4r3fn0c * 4c 30r 4c « « « 4c « « 3ftc 3^ 4c 3ft 4c t ♦ « 3^ 4 c 3» ♦ ♦ 4c ♦ ♦ 3^ « «4C 4 c3|c « 3^ « « 3^ * 4c « « 4c 


RrAO( «V) 

( W( n t 1 = 1 t NN) 

REWIND 

A 

on Q2 J = 

1 » 3 

U ( J ) 

O.DO 

nn 92 I = 

1 , NN 

KKK “ 

(J-l)^NN ♦ I 

1? ( J ) 

U( J ) - A(KKK)«W( I ) 

CONTINUE 

CALL DGMPRD(NI,U,DX,3,3,1) 

nc 9b 1 = 

1 f 3 

JK = 

( I-1)«NN + 1 

JM 

JK + NN - 1 

RF4D(3> 

(A(J), J= JK , JM) 

P.FWIND 

3 


4oft3ft4 COMPUTING VARIANCE OF UNIT WEIGHT • M02 • 


DO 9/S I ~ 

W( n 

no 9f J = 

K 

W( 1) 

96- Cr?NTINUF 

RFA0(4) 
REWIND 
DO 97 K = 
W(K) = 

^7 CONTINUE 
MMM = 

nr; 98 k = 


1 » NN 
O.DO 
1 » 3 

(J-1)*NN + I 

W( I ) - A(K)*DXU) 


(A( 1) , 
A 


I , NN 
W(K) - 


U 1 


A(K) 


NN) 


0 

1 


NN 
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A(K) 

= 

O.DO 

LI 

= 

( (K-l)/2)*2 ♦ K 

L2 

= 

LI * 2 

DC; 98 L 

- 

LI t L2 , 2 

MMM 

= 

HL-l)/2) ♦ X 

AIK) 

= 

A(K) * P(L)4W(MMM) 

»FAD(9) 


(Wm , I = 1 , NN) 

REWIND 


4 

VPV 

- 

0.00 

DC: 9P K 

= 

1 * NN 

VPV 

= 

VPV - A(K)*W(K) 

M02 

=■ 

VFV/INN - 3) 


C 

C 

c 

C CCMPUTING VARIANCE^ COVARIANCE MATRIX • VAR 

C 

c 


100 


105 


DO ICO 1=^ 
DO ICO J= 
VAR( I T J>- 
CC-NiTlNUE 
nn 10*5 1 = 
DX (I ) 


1 , 3 
1 » 3 

Ma2*RH0S*RH0S^NI ( I f J I 
1 ♦ 3 

DX( I )*RHOS 


t 


♦ COMPUTING COEFFICIENTS OF CO-RELATIONS FOP PARAMETERS 


on 110 1= 

IF 

JJ = 

on 106 j= 

Mn,j) = 

106 NI ( j, n = 

1G7 NKIfl) = 
110 CONTI NUF 


1 ♦ 3 

(I. EC. 3) GO TO 107 
I + 1 
JJ f 3 

VAR(I,J»/(DSORT{VAR( I,n » ♦DSCRT( VAR ( J, J HI 
NI(I,J) 

i.no 


4 ^ ^ ^ Ac :0c « 4c 4> ♦ ♦ « 4c* 

. * 
♦ 

**A* WRITING OF FINAL SOLUTION VECTOR AND VARIANCE-COVARIANCE MATRIX ♦ 

* 

♦ 

4* 4 A *****4 **************** *************************************************** 


WP ITF(-6,60?5 ) 

6025 format (• 1 • t///) 

WPlTe(f ,6028) (N»ME1{ I) , 1 = 1»3) f (NA MP2 < I) 1 1 = 1 , 3 1 

602? FPPKAH' • ,5Xt3A4, '-TO-' f3AAf/f 

?6X , •#*»:♦**♦♦♦*♦*♦*♦*♦*♦*♦♦«**♦♦* * / f / n 
WPITPt6,603C) 

6C30 FORMATC • ,30X , • SCILUTI ON FOR "3** ROTATION PARAMETERS',/, 

231X, » ',/, 

32'^-X, • (FROM DIRECTION COSINFS — UNITS SECONDS OF ARC)',./) 
r-n Tn n 12 , 114 , 116 ), kpr 
112 WR1TE(6,6031) 
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(‘•031 FORMATC • f37 X, musing VARIANCFS ONLY)',//) 

GO TO 120 
llA WR ITF (6,6032) 

6032 FORMATC ',21X, 

2* (USING 3X3 BANDED DIAGONAL VARIANCE-COVARIANCE MATRIX)',//) 
GO TC 120 
116 WRITE(6,6033) 

6033 FPPMATC ' ,29X, '(USING FULL VARIANCE-COVARIANCE MATRIX)',//) 
120 WR1TF( 6,6035) 

6035 FDRMATC ' , 20X , 'OMEGA • , 19X , ' PS I ' , 20X, 'EPSI LON* , //) 
WP1TE(6,6040) (DX( I ), 1=1,3) 

6040 FORMAT! ' ', 5X , 3D24.7 , // ) 

WP1TR(6,6045) 

6043 FORMAT! ' ' ,32X , » VAR I ANCE - COVARIANCE MATRIX',/, 

233X, ' ',/) 

WRITE(6,604R) M02 

6041'. FORMAT! ' ' , 1 7X , ' M02=' , F6 . 2 , // ) 

WPITE(6,6050) ((VAR(1,J), J=l,3), 1=1,3) 

6050 FORMAT!' ', 3X ,3025.8 ,// ( 4X , 3025. 8 , /) ) 

WPITF!6,6075) 

6075 FORMAT!' • ,33X , 'COEFF I C 1 FNT OF CORR FLAT ION* , /, 

234X, • ',//) 

WRITE (6,6085) ( ! NI ( I , J ) , J= 1 , 3 ) , 1=1 , 3) 

6085 FORMAT!' ', 3X , 302 5 .6 , // ! 4X , 3025. 8 ,/ ) ) 

IF (KCODE(ll) .FO. 3) GO TO 150 

WPITF (6,7000) 

7000 FPCMAT !'-' ,//,34X, 'SOLUTION FOR SCALE FACTOR',/, 

234X, ' ' ,/, 

335X,'!FROM CHORD COMPARISON)',//) 

WRITE (6,7004) 

7004 FORMAT!' »,20X, 'SCALE FACTOR ', 27X ,' VA R I ANCE' ,/, 

22? X, • (10.0+5) ' ,29X, ' llC.D+11) ',//) 

WPITE (6,7010) S , VSF 

7010 FORMAT!' ' ,20X ,F8 . 2,30X, F7.2 , // ) 

150 KCOOF(ll) = 4 

RETURN 
END 
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SUBROUTINE TFORM (NO, NO) 

1^^PL1CIT REAL ♦ 8 (A-H, G-Z ) 

REAL ♦ 8 MI ,KK,KL,NI,M02 

n I y,E NS IDN XYZ (35,3) ,UVW( 35,3) , S IGMAX( 7,7) , NAMEl ( 3) , 

2A( 3600) ,W( 1200),VAR(7,7) ,DX(7),NI (A9),NSTA(35) ,NAME2(3) , 
:ANG(A) ,U(7),LT(7),KC0DEa5),CNT(7,4),TT(7,4) ,CN(4, 7) , ZP ( 4,4) , 
AM.T(7) ,KSTA(35) ,VP(7,7) ,X0 (7) , KL (150 ) , KK ( 150) ,MI ( 2400 ) ,ROT( 4,4) 
COMMON /WEIGHT/ MI 

COMMON /CODE/ KCODE 

COMMON /ANGLE/ ANG,ROT 

COMMON /INAME/ NAMF1,NAME2 

COMMON /CRNT/ VPV , OX , S02 , XD , S IGM AX 

COMMON NSTA,KSTA,NN,NM,UVW,XY2,A,W,KPR,KPARM 

Pil= 3.141592653589793DC 
RHC = 180.D0/P1I 
RHOS = RH0+3600.D0 
I PAR A = KC0DE(2) 

IC = KCOOE(ll) 

KPUNT = 1 
DO 5 I = 1 , 4 

no 5 J = 1,7 

CN(I,J) = 0.00 

TT(J,I) = 0.00 

5 CNT(J,I) = 0.00 

nn 10 I = 1,4 

DC 10 J = 1,4 

10 zP(i,j) = o.no 


SETTING UP MATRIX *A» — COMMON TO ALL SOLUTION 

C 

C 

c 

NNN = 6* NO 
NN2 =: NO^IPARA 

13 I = 1 , NNZ 
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AM ) = O.DO 

CHNTIMUE 

no 15 1 = It NO 

KKK = (3*l-2) 

LLL = KK.K4NQ+1 

MMM LLL+NO+1 

A^KKK) = 1 .DO 

A(LLL) = l.PO 

A(MMM) = 1.00 


SETTING UP MATRIX »W» WHICH IS COMMON TO ALL SOLUTION 


W(KKK) = (UVWdtD-XYZdtl)) 
W(KKK+n= (UVWdt2)-XY7(It?)) 
W(KKK>2)- (UVWd ,3)-XY2( 1,3) ) 

15 CONTINUE 

IP (KCQDE(2) .NF. 3) GO TO 50 


SOLUTION FOR 3 TRANSLATION PARAMETERS 


N = 3 

IC^SF = 1 
CO TO 81 


SOLUTION FOR 3 TRANSLATION AND 1 SCALE PARAMETERS 

C 

C 

SO N - ^ 

DC 60 I = 1, NO 
KKK = 3*<N0+l)-2 

A(KKK) = UVW( 1,1) 

A(KKK+1)= UVW( 1,2) 

A(KKK + 2)= UVWd,3) 

60 CCNTINUF 

IF (KC0DF(2) .NE. 4) GO TO 70 
ICASF 2 

GC TO 61 

C 
C 

c 

0 SOLUTION FOR 3 TRANSLATION, 1 SCALF AND 3 ROTATION PARAMETERS 

C. 

C 

r 

70 N = 7 

ICASF = 3 

on 60 1 = 1, NO 

KKK = 4«N0-*-(3*I-2) 
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LLL = 

KKK + NO 

MMM = 

LLL 4 NO ♦ 1 

A(KKK.) 

= UVW( I ,2) 

A(KKK4»1) 

=--UVW( 1,1) 

A(LLL) 

=-UVW(I,3) 

A(LLL+2) 

= uvwn,i) 

A (MMM) 

= UVW(I,3) 

A'(MMM4-1) 

=-UVW(I,2) 

CONTINUE 
on I 

= 1 , N 

KKK 

= (I-D+NQ+l 

LLL 

= KKK4-NQ-1 

WRITF (3) 

(AU), J=KKK,LLL> 

CONTINUe 

RFWIND 

3 

WPITB( A) 

(W(I), 1=1, NO) 

RFWIND 

A 


F^'?M1^^G NDOMAL eOUATICNS — MAT*?ICFS 'N* AND 'U* 


« 4c « 4c 4 « 4c 4c « « 4^ 4c «* « 4C4 4c 4c 4c « 4C« ^ 4c :^4c « «4c « « 4c4c ♦♦ 4 c4m^ ♦ « #4C 4c4c4i ««4c 4c« 

100 CALL SFTUP (^0*NO,IPAP.A) 

DC IIB I = 1 , M 
RFA0(?) (W(J)» J=1,N0) 

K1 = (I-n*N0+l 

K? = Kl+NQ-1 

MMM = 0 

DO 116 K= K1 , K2 

A(K) = 0,00 

LI = (( (K-Kl)/3)«3)-U 

L?. ^ LI + 2 

DC 116 L = LI , L2 

MMM = MMM + 1 ■ 

116 A(K) = A(K) 4- W(L)*MI(MMM) 
lie CCNTINUF 

RFWIND 3 

DO 120 I = 1 , N 
'>rAD(3) (W(L), L= 1,N0) 

JK = (i-n*w + i 

JL = JK+M-l 

DO 11^ J = JK,JL 
M ( J) = O.DO 
DO 119 K = 1 , NO 
II = (J-JK)*NQ -*■ K 

no NKJ) = NKJ) + A(lI)4cW(K) 

12C CONTI NUF 

RFWIND 3 

no 121 I = 1 , N 
no 1?1 J = 1 , N 
K = ( I-1)«N +J 

121 SIGMAX(I.J) = NHK) 

°FAD(6) (W(I)» I= ItNO) 

RFWIND A 
on IZ2 J=1 t N 
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UU) := O.DO 
DO 122 I=l , NO 
KKK = ( J-1)*NQ+I 

U( J) = U(J) « A(KKK)*W( I) 
12? CONTINUF 


* 

CQ^^PUTING SOLUTION VECTOR »DX» FOR TRANSFORMATION PARAMETERS ♦ 

♦ 

4c :4c; 4c ^ ifc 4c 4c4c :4c A 3^ 4c « 4c(c:4r 3«t 4c « 4c# 4c 4r 4c « 4 4c ^ 4i;0c * 4(# 4c# 4( « 41 * 4i ^ 

CALL DMlNVtNl tNtDTtLT »MTI 
CALL DARRAYI 1,N,N,7,7,NI »VR ) 

CALL OGMPRn(NltU,DX,N,Nf 1 ) 
on 123 I = 1 , N 
JK = (I-l)*NO 4 1 

JM = JK 4- NO -X 

123 '^FA0(3) ( A{ J) tJ^JKrJM) 

P^'WIND 3 


###* COMPUTING VARIANCE OF UNIT WEIGHT •M02* 


1?3 


126 


no 123 

I 

= 1 f NO 

W( I ) 

n 

O.DO 

no 123 

J 

= 1 f N 

K7X 


= (J-1)*N0+I 

W( I ) 

= 

wm -A(KZX)+DX(J) 

CONT 1 NUF 


READ! A ) 


I A I n , I = 1 , NO ) 

Pf'WiND 


4 

no 126 

K 

= It NO 

W(K) 

= 

WIK)-A(K) 

CHNTINUE 


MMM 

= 

0 

on 12B 

K 

= 1 f NQ 

AIK) 

= 

O.DO 

LI 

( (K-l ) /3)=f6 + K 

L2 


LI 6 

DC 12E 

L 

= LltL2,3 

MMM 


= I(L-l) /3) +1 

AIK) 


= A(K) ♦MI(L)*W(MMM) 

CALL 


RFSIDU IN0,NNN) 

PEAHI A) 


(W( I ) »I= 1,N0) 

REWIND 


4 

VPV 

= 

O.DO 

no 130 

K 

= ItNO 

VPV = 


VPV ^ AIK)#WIK) 

MO?. = 


VPV/INO-N) 


12R 


130 


C 

c 

r 444# COMPUTING VAR IANCF-COVAP I ANCE MATRIX 'VAR* 

C 

C 

on 13? I - 1 , N 
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no 132 J = 1, N 
VAP.< I , J) = M02*VR( I, J) 
132 CONTINUE 


IF (KCCDE(2) .EO. 

3) 

GO 

TO 

140 

DX(4) = DX(4) ■•= 10 

.05 




IF (KCDDF(2) .EO. 

4) 

GO 

TO 

140 


DO 135 I = 5 , 7 
DX( n = DX( I) * RHOS 
135 CONTINUE 

**** COMPUTING COEFFICIENTS OF CORRELATIONS FOR PARAMETERS 


L4C 

DO 145 I = 

1 fN 


IF(I .EO.N) 

GO TO 144 


JJ = I ♦ 

1 


DO 142 J = 

JJ » N 


VR(1,J) = 

VAR( I,J)/(DSORT(VAR( I , I ) ) ♦DSORT ( VAR ( J# J) )) 

142 

VR(J,I) = 

VR(ItJ) 

144 

VR(I, 1 ) = 

l.DO 

145 

CONTINUE 


200 

WPITE((^»250) 

250 

FORMAT! ' 1' 

t//) 


WRITF(6, 3001 (NAMFII I ) ,I=lf 3) t (NAMF2( n ,1=1,3) 
300 FORMAT) • • ,5X,3AA, '-TO-' ,3A4, /, 

26 X 

60 TO (500,600,700) , ICASE 


:;(4:****’i^*:**************4^**********************************************«*****> 


*♦«* WRITING OF FINAL SOLUTION VECTOR AND VARIANCE-COVARIANCE MATRIX 




500 WRITE (6,6025) 

6025 FORMAT (• *,/////) 

WRITE(6,6030) 

6030 FORMATC-' ,21X, ‘SOLUTION FOR 3 TRANSLATION PARAMETERS',/, 
232X, '(UNITS - METERS)',///) 

GO TO (512,514,516) , KPR 
512 WP1TE(6,6032) 

6032 FORMATC ' ,29X ,' (USING VAR lANCES ONLY) ',// ) 

Go TO 520 
514 WRITE(6,6034) 

6034 FORMATC ',15X, 

?• (USING 3X3 BANDED DIAGONAL VARIANCE-COVARIANCE MATRIX) ',//! 
GO TO 520 
516 WR ITE(6,6036) 
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f.036 FGRMfrTC • ,22X, • (USING FULL VARI ANCF-COVAR 1 ANCE MATRIX) ',//) 

520 WRITE(6,6038) 

6036 format ( fl6X,*nx*,20X, 'DY* ,22Xf 'DZ* ,//) 

WRITF(6,6040) (DX(I), 1=1,3) 

6040 FORMAT) ' IX, 3D23. 8, //////) 

WRITE (6, 6045) 

6045 FORMATC ' ,26X , • VAR 1 ANCE - COVARIANCE MATRIX*,//) 

WPITF( 6,6048) M02 

604G FORMAT ( , 14X , • M02= • , F6 . 2 ,// ) 

WRITF(6,6050) ((VAP(I,J), J=l,3), 1=1,3) 

6050 FORMATC 1 X ,3023.8 ,// ( 2X , 3023 . 8 , / ) ) 

WRIT'=-(6,6075) 

6075 FORMATC-' ,//,27X, 'COEFFICIENTS OF CORRELATION*,////) 

WRITF(6,6085) CVR( I,J) , J=1,N), 1=1, N) 

6085 FORMATC IX ,3023.8 ,// ( 2X , 3023. 8 , / ) ) 

GO TO 1000 
600 WRITE(6,6500) 

6500 FORMAT! * *,/////) 

WR1TE(6,6510) 

6510 FORMAT!'-* ,17X, 'SOLUTION FOR 3 TRANSLATION AND 1 SCALE PARAMETERS* 
2,/,34X,'!UNITS - METERS)*,///) 

GO TO !612,614,616) , KPR 
612 WRITE (6,6512) 

6512 FORMAT!' • ,29X ,'( USING VARIANCES ONLY)',//) 

GO TO 620 
614 W» 1TE(6,6514) 

6514 FORMATC *,19X, 

2' (USING 3X3 BANDED DIAGONAL VAR I ANC F-COVAR 1 ANCE MATRIX)*,//) 

GO TO 620 
616 WP1TF(6,6516) 

6516 FORMAT!' ' ,22X, • (USING FULL VARIANCE-COVARIANCE MATRIX)*,//) 

620 WPITE!6,6520) 

6520 FDPMATC ', 6X , * DX * ,22X, » DY • , 23X , • 02 * , 22X , * DL * ,//) 
WRITe!6,6550)!DX!I), 1=1,4) 

6550 FORMATC ',015.8,3024.8//////) 

WR1TE!6,6600) 

6600 FORMATC • ,26X, • VARIANCE - COVARIANCF MATRIX*,//) 

WRITF(6,6625) M02 

6625 FORMATC 8X , • M02=* , F6 . 2 , // ) 

WRITE!6,6650) !!VAR(I,J), J=l,4), 1=1,4) 

6650 FORMATC *, IX ,4020.8 ,// ! 2X ,4020 .8 , / ) ) 

WRITE(6,6675) 

6675 FORMAT! •-• ,//,27X, 'COEFFICIENTS OF CORRELATION*,////) 

WR1TE!6,6685) ! !VR| I,J) , J=1,N), 1=1, N) 

6685 FORMAT!' IX ,4020.8 ,// ! 2X ,4020.8 ,/) ) 

GO TO 1000 

70C GO TO ! 710,705) , KOUNT 
705 IF (KPARM .EO. 4 ) GO TO 708 
WPITE (6,7002) 

7002 FORMATC • ,28X , • ROTAT I ON PARAMETERS CONSTRAINED*,/, 

GO TO 710 
70S WR1TF!6,7005) 

7005 FORMAT!* *,20X, 'SCALE FACTOR AND ROTATION PARAMETERS CONSTRAINED*, 
710 WP1TE!6,7010) 

7010 FORMAT!' ', 13X ,* SOLUTION FOR 3 TRANSLATION, 1 SCALE AND 3 ROTATION 
2 PARAMETFRS* , /,14X, ' 
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2,T58,' •,/) 

GP TO (712,714,716), KPR 
71? WRITE(6,7012) 

7012 FORMATC • ,34X, • (USING VARIANCES ONLY)*,//) 

GO TO 720 
714 WRITE (6,7014) 

7014 FORMAT ( • • ,16X, 

2* (USING 3X3 BANDED DIAGONAL VARIANCE-COVARIANCE MATRIX)*,//) 

GO TO 720 
716 WRITE(6,7016) 

7016 FORMATC ' ,24X, * (USING FULL VARIANCE-COVARIANCE MATRIX)*,//) 

720 WRITE (6,7020) 

7C20 FORMATC ',16X,*DX*, 6X,*DY», 6X,»DZ', 7X,'DL*, 5X,*0MEGA*, 

2T S9,*PSC, 4X,*EPSIL0N*,/, 

3 IS X, 'METERS* , 2X,* METERS*, 2X , • ME TFRS • , IX, * ( 10 .0+5 ) * , IX , * SECONDS* , 

4TE7, • SECONDS • , 2X ,* SECONDS *, /) 

WRITE- (6, 7030) DX 

7030 FORMATC • , 1 2X , F7 . 2 , 2F 8. 2 , F 8 . 2 , F9 . 2 , T55 , F8 . 2 ,F9 . 2 , // ) 

WRITE(6,7040) 

704G FORMATCO* ,28X, 'VARIANCE - COVARIANCE MATRIX*,//) 

WR ITE(6,7045) M02 

7045 FORMAT!' ' , lOX , * M02= * , F6 . 2 , / ) 

WP 1TE(6,7050) ((VAR(I,J), J=l,7), I = 1,7) 

7050 FORMATC ' , 2X , 7D1 1 .3 , / / ( 3X ,7D1 1 . 3 , / ) ) 

WP)TF(6,7075) 

7075 FORMAT!' *,/, 29X , 'COE FF IC I ENTS OF CORRELATION*,//) 

WRITF(6,7085) ( (VR( I,J) , J=1,N), 1 = 1, N) 

7085 FORMAT!' ', 2X ,7D1 1 . 3, // ( 3X , 7D11 .3, / ) ) 

IF(IC.EQ.O) GO TO 1000 
WRITE (6,7090) 

7090 FORMAT (* 1 ',////, 36X ,' RES I DUALS V',/,36X,* *,///, 

212X, 'FIRST SYSTEM' ,33X, 'SECOND SYSTEM*,/) 

KSM = NNN ♦ 1 
KMR = NNN - 1 +KSM 

WRITE (6,80CC) (A(I), I = KSM, KMR) 

8C00 FORMATC * ,4X , 3F8 .3, 22X, 3 F8 .3 , / ( 5X, 3F8 .3 , 22X , 3F8 .3 ) ) 

IF (KC0DE(3) .EO. 0) GO TO 1000 

C 

(; ,>^9^ A ^4:**^*’^*^’^**'t‘**'**^4’l‘*l^ ********’************************************«****** 


c ♦ 

C * 

c ♦ 

c OBTAINING CONSTRAINED SOLUTION FOR ROTATION PARAMETERS * 

C * 

C ♦ 

C * 


c 

c 

C^LL CSTRNT(N,NQ, IC,U,CN,CNT,TT*2P) 
KCCCFO) = 0 

on 725 1=1,7 

nxn ) = xD(T ) 

DC 725 J = 1 , 7 
V&R( 1 , J) = SIGMAXC I,J) 

725 CPMTIMUE 

on 750 I = 1,N 
IF(I.eo.N) GO TO 740 
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JJ =1+1 
DO J = Jj , N 

VR n , J) = VAR( I, J)/(DSQRT( VAR( I , I )) SORT ( VAR ( J» J) ) ) 
735 VR(J,U = VR(ltJ) 

740 VR(I,I) = l.DO 
730 CONTINUE 
KOUNT = 2 
MG2 = SOZ 

GO TO 200 
icoo return 

END 
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c 





c 





c 





c 

♦♦♦♦♦♦♦♦«<<«««♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦«««««♦♦♦♦♦♦♦♦ > 

c 




*«*«**#« 

c 





c 



SUBROUTINE SCALE 


c 





c 





c 





c 

♦ ♦♦♦♦♦♦♦«««««♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦*♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦»««««<<♦♦♦♦♦♦♦♦ 1 

c 





c 




«**«#*«« 

c 



FINDING WEIGHTED MEAN AND VARIANCE FOR 


c 



•SCALE FACTOR* BY COMPARISON OF CHORDS IN 


c 



THE TWO SYSTEMS BY CALLING SUBROUTINE »SCALE*. 


c 





c 





c 





c 

♦ ♦♦♦♦*♦♦«««««♦♦♦♦♦»♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦#♦♦♦♦♦«««««♦♦♦♦♦♦♦♦ 1 

c 





c 




*«*«««*« 

c 





c 






SUBROUTINE 

SCALE (NOtNtStVSFfWT) 



IMPLICIT 


REAL • 8 (A-H i 0-Z) 



DIMENSION 

P(12tl2).H(12l,PF(6,6),PS(6t6), 



2 


HI(12)fDL(600) vVI(600) >WI(600) 



COMMON 


/SFAC/ SW,SF,DUtDV,DW,DX|DY,DZ,Rl,R2tPF,PS 



RR 

s 

R1 * R2 



RT 

s 

R2 /IR1**3» 


c 

<c4t««#4r^4c 




c 





c 





c 


SETTING UP OF VARIANCES FOR EACH CHORD THROUGH ERROR 

PROPOGATION 

c 





c 





c 






H(l) 

C 

DU ♦ RT 



H(2) 

= 

DV ♦ RT 



H(3) 

= 

DW ♦ RT 



H(7) 

s 

-OX/RR 



H(8) 

= 

-OY/RR 



H(9) 

= 

-OZ/RR 



DO 10 I 

= 

1 » 3 



H( 1 + 3) 

= 

-HU) 



H( I + 9) 

= 

-H(I ♦ 6) 



10 CONTINUE 





DO 15 I 

= 

1 f 12 



DC 15 J 

= 

1 t 12 



15 P(I,J) 

= 

0*00 



DO 20 I 

s 

1 f 6 ; 



L 

= 

I ♦ 6 



DO 20 J 

= 

1 t 6 



M 

= 

J ♦ 6 



P(It J) 

= 

PP(I,J) 



20 P(L,M) 


PS(I,J) 



CALL 


DGMPRD (H,P,HI,l,12fl2) 
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CALL 

DGMPRO (Hl,H,WS,ltl?,l) 

c 



c 



c 


FINDING WFIGHTED MEAN FOR SCALE FACTOR 

c 



c 




WS 

= l.DO/WS 


WI (N) 

= WS 


SFl 

= R?/RX - 1 


DL(N) 

= SFI 


SF 

= SF ♦ SFI ♦WS 


sw 

= SW •». WS 


s 

= SF/SW 

c 



c 

j^:Of « + 


c 

:<c5}r4t5jc«:»£«^* 


c 


FINDING VARIANCE FOR THE WFIGHTED MEAN 

c 



c 



c 




IF 

(N .NF. NO) GO TO 500 


PVV 

= O.DO 


DP 50 K 

- 1 t NO 


VHK) 

= ( (S-OUK) )**?)*WI(K) 


^0 PV\/ 

= PVV + VI (K) 


VSF 

= Pvv/(sw*(Nn-i) ) 


5 

= S ♦ 10. D5 


WT 

= l.DO / VSF 

r 

.00 R*^TURNi 
E=NiO 
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« « 4 « 4c 34 c ^ ^ « ♦ # # ♦ ♦ ♦ « ♦ « 4c# « 4 # « « 

♦ 

• CSTRNT » 0 

♦ 

4c 34 c 3 ^ #3«c 3fc 4^ :4c 4(4^4c 4c# 3 ^ ###3^#«3Cc#34^«3g(#34c#3|C##3(C 3 ^ 3«( 30c 3|13 Cc#3»##4C##«#### ##############«###«###### 

♦ 
♦ 
# 
# 
♦ 
4i 
♦ 
♦ 
♦ 

########################################################«#######♦#####«#####* 
######## #4(###### # 

#4i###### ####♦♦♦♦ 4t 

SUBROUT INF C STRNT ( N»NN » I C • WS , CN t CNT , TT » ZP ) 

IMPLICIT PEAL ♦ 6 (A-H , D-Z) 

REAL ♦ 8 M02,KC 

DIMENSION XD(7) *WS(7) ,WX(7) •KU^) , 

?WC(4) ,LM(7),MM(7) fPZ (4,4) ,CNT(7, IC) ,GG(7,7), 
:SIGMAX(7,7),nX(7)»TK(7),TT(7,IC),CN(IC,7),ZP(ICf IC) 


34.#44r SUBROUTINE SOLVES FOR CONSTRAINFD CASE IN RESPECT OF *3» 

4t##4i ROTATION PARAMETERS. CONSTRAINTS ARE CODED FOR •ALL* THE 
#### PARAMETERS — BLANKS CARDS ARE NEEDED FOR NON-CONSTRAINTS. 

INPUT CONSTRAINTS FOR ROTATION PARAMETERS ARE IN SECONDS OF ARC. 


23 


ccmcu 

COMMON 
PI I 
RHD 
PROS 
WC(1) 
on 23 
WC ( I) 
CONTINUE 


I = 


/CPNT/ VPV,DX,S02,XD,SIGMAX 
/angle/ WC,PZ 

3.141592653589793D0 
= IBO.DO/PIl 
= RHO ♦ 3600.no 
WC(1) / 10.05 
2 , IC 
WC( I)/ RHOS 


#### SETUP CONSTRAINTS MATRIX •CN* REQUIRED FOR SOLUTION 



DP 25 I 


1 

. IC 



DO 25 J 

= 

1 

» N 



CNN , J) 

= 

0. 

»D0 



TT( J , 1 ) 

= 

0, 

.DO 



CONTINUE 






IF ( IC . 

EO 

• 

4) 

GO 


no 50 i 



1 f 

IC 


DC 50 J 

= 


1 » 

IC 

EC 

?P(1 ,J) 

‘ = 


PZ(I 

♦ 1 


CC TO 

200 



ICG 

on 150 I 



1 » 

IC 


or 150 J 

= 


1 f 

IC 

150 

ZP(1, J) 

=: 


PZd 

t J 

2C0 

no 300 I 

=■ 

1 

f 4 



J 

= 

I 

+ 3 



CNN , J) 

r- 

1. 

.0 


3 or 

CCNTINUe 
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c 



C ^CLVe FOR EFFECTS OF CONSTRAINTS ON THF SOLUTION VECTOR 'OX* 

C neTAINFO FRO?^ NON-CONSTRAINT SOLUTION 



C 

no 520 1= 1 , 1C 
00 520 J= 1,N 
CNT(J,I)= CN(I,J) 

520 CONTINUE 

CALL MTPY(CNT,ZP,NtIC, IC,TT) 

CALL MTPY(TT,CN,N,IC*N»&G) 
on 522 I = 1 , N 
DO 522 J = 1 ♦ N 

522 GG(I,J) = SIGMAXCItJ) ^ r,G(I,J) 

CALL nMlNV(GOtN»DTT,LM,MM) 

C£LL MTPY(TT,WC,N, IC,1,WX) 

00 525 1= 1 t N 

wsm = (wsm - wxnn 

525 CONTINUE 

CALL PTPYIGGf WS,N,N,1,X0) 

C 

C COMPUTE NEW VARIANCE OF UNIT WFIGHT AND 

C ^•F;w VARIANCE - COVARIANCE MATRIX 

C 

CALL MTPY(CN»XD,IC»N,1,KC) 

00 535 I = 1 , IC 
535 KC n ) = -KC( I )-WC( I ) 

CALL MTPY(P2#KC»IC,IC»lfDX) 

SUM = 0.0 

DO 5 AG 1= 1, 1C 

SUM = SUM + DX ( I ) ♦ wc m 

5A0 CONTINUE 

?VV = VPV - SUM 

SC2 = PVV/(NN-N+1C ) 

no 550 1= 1 , N 

DO 550 J= 1 , N 

SIGMAXdfJ) = S02*GG(I,J) 

550 CONTINUE 

Xn(4) = xn(4) ♦ 10. D5 

DC 560 1=5,7 

X^(I) = X0( 1) ♦ PROS 

•^60 CONTINUE 
1000 RETURN 
END 
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o n r< 


* 

SETUP * 

* 

* 

SETUP MATRIX B»TP — Ml — AFTEP READING VARIANCE ♦ 

COVARIANCE MATRIX FOR EACH POINT SEPARATELY AND ♦ 

TH'^N STORING THF ELEMENTS SO FORMED IN THE PROPER PLACE IN >MI». ♦ 

♦ 


SUBROUTINE SETUP tNO,NN,lPARM 
IMPLICIT REAL ♦8<A-H,0-Z) 

REAL E MINK 

OIMENSION R(3»6) ,RT{6,3) ,PI(6,6) ,PK(3,3) t 
2XM(3 t? ) ,XK (3,6) ,MINK( 2AOO) ,LM (3 ) , MM ( 3 ) , KCODE ( 15 ) 
COMMON /RES/ RT 

CrMMPN /WFIGHT/ MINK 

CC:MM(}N /code/ KCOOE 


SETTING UP MATRIX 'B* WHICH WILL BE SAME FOR ALL SOLUTION 

C 

C. 

NV = NN^3 

OG 8 1 = 1 , NV 

MINKd) = 0.00 
8 CnNTINUF 

DO 1C 1 - 1 , 3 

DO 10 J = I , 6 

F( 1, J) =0.00 
XK ( I , J ) =0.00 

10 CONTINUE 

B(l,l) = -l.DO 

P.(2,2) = *1.00 

B(5,3) = -l.DO 

8(1, A) = l.DO 

P(?,5) = l.DO 

8(3,^) = l.DO 

DP 12 1 = 1 , 3 

DP 1? J - 1 , 6 

RT(J,1) = B(I,J) 

12 CONTINUE 

OP 15 1=1,6 
nn IF J = 1 , 6 
t>I ( I ,J ) = O.DO 
15 CPNTINUE 

no ?C I X 1 , 3 

nr 20 J = 1 , 3 

XM( I , J ) = O.DO 

20 CP'NTINUF 

KMS = 0 

IP (KC0DF(14) .EO. 1) GO TO 65 
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IF (KCODEie) .EQ. 1) GO TO 54 

00 40 L = 1 t NO 

DC 39 J = 1 t 3 

READ(5,38) (PI(J,K), K=l,3) 

38 FORMAT(3F5.2) 

39 WRITE (2) (PI(J,K) f K = I , 3) 

40 CONTINUE 

DO 52 M = 1 , NO 

DO 44 J = 4 , 6 

RFAD(5,42) (PI(J,K), K= 4, 6) 

42 FORMAT(3F5.2) 

44 CONTINUE 

DO 45 I = X , 3 

DO 45 J = 1 , 3 

45 PK(I,J) = PH 1 + 3, J-t-3) 

DO 50 I = X , 3 

50 WRITE(X) (PK(I,J), J=Xf3) 

52 CONTINUE 
REWIND X 
REWIND 2 
GO TO 65 

54 DO 58 M = X » NO 

READ (5,55) (PHI,I) , 1= X,3) 

55 FORMAT (3FX0.2) 

DO 56 I = X , 3 

56 WRITE (2) (PI( I,J) . J = X,3) 

58 CONTINUE 

DO 64 M = X , NO 

READ (5,55) (PHI, I) , 1= 4,6) 

DO 60 J = X , 3 

DO 60 K = X , 3 

60 PK(J,K) = PI(J+3,K+3) 

DO 62 1 = X , 3 

62 WRITE(X) (PK(I,J), J = X,3) 

64 CONTINUE 
REWIND X 
REWIND 2 

65 DO 100 I = X , NO 

KKS = ( I-l)*9 ♦ X 


**** READ IN VARIANCE - COVARIANCE MATRIX AS BLOCK DIAGONALS 
*44* OF (6,6) MATRICES FOR EACH POINT USED IN TRANSFORMATION. 
*444 MATRIX 'PI* IS BUILT UP POINTWISE - FIRST (3,3) BLOCK 

REFERS TO SECOND COORDINATE SYSTEM AND SECOND (3,3) BLOCK 
4*4* THEN CORRESPONDS TO FIRST COORDINATE SYSTEM, 


DO 70 J = X , 3 

READ(2) (P1(J,K), K= X,3) 

70 CONTINUE 

DO 74 L = 4 , 6 I 

READ( 1 ) (PI (L,M) , M=4,6) 

74 CONTINUE 

CALL MTPY(B,PI,3,6,6,XK) 

CALL MTPY(XK,BT,3,6,3,XM) 

CALL DMINV(XM,3,DET,LM,MM) 

MINK(KMS ) = XM(X,X) 

- 60 - . 



100 


MINK(KMS+1 ) 

MIWK(KMS^2) 

MINK(K«S+3) 

MINK(KMS+4I 

MINK(KHS>5) 

MINk(KMS^6) 

MINK(K«S^7) 

f^INK(KMS4-8) 

CQN7 I NUE 

REWIND 1 

REWIND 2 

PFTUPN 

END 


XMI3,1) 

XHUt2I 

XM(2f2) 

XM(3,2) 

XMIltS) 

XHI2,3) 

XM|3f3) 
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«««*«**« « 

]|[ 

4;0>!0L,9c4r«i^* ««*«««*♦ ♦ 

**:fr*«^«* SUBROUTINE »RES10UF» ♦♦♦♦♦♦♦♦ * 

«*«**««« * 

7^*H^:i^^i^it^i^<<««««**^^*^**^*^****^*******^********^*^*****«««««^***^^** ♦ 

^ 4 c ^«; 0 c ^« 4 c **♦♦«#«* « 

♦♦♦«♦♦«« this subroutine computes residuals for each ♦♦♦♦♦♦♦♦ * 

SYSTEM COORDINATES (USED AS OBSERVATIONS I • ♦♦♦♦♦♦♦♦ « 

♦ 4i 


:^:4c«4c:^A:^4c 
T. ♦**«3fr*«* 

SUPROUTINE 

IMPLICIT 

DIMPNSION 

? 

COMMON 
COMMON 
COMMON 
DP 5 I = 
DO S J =: 
^ PKIfJ) = 
DC 25 1 = 

JJ 

DC 10 J = 
L 


«««*#««« « 

*«««*«*« 4( 


RESrOU (NO » NNN) 

REAL ♦ 8 (A-H t 0-Z ) 

NSTA(35),KSTA(35),UVH(35,3),MM(6I,BT(6.3),LM(6), 
XY2(35,3) ,A ( 3600 ) ,H( 1200) tPI (6 »6 ) tKCOOE ( 15 ) ,8S( 6 »3 ) 
/RES/ BT 
/CODE/ KCODE 

NSTA»KSTA,NN,NM»UVW,XYZ*A,W 
1 ♦ 6 
1 f 6 
0*00 
1 , NO 

NNN ♦ C I-l)«6 
1 f 3 
J ♦ 3 




READING VARIANCE - COVARIANCE MATRICES FOR 
THE FIRST SYSTEM ~ POINTWISE ~ AS (3X3) . 


prAD 


(1) (PI(L*M)t M = 4f6) 




READING VARIANCE -- COVARIANCE MATRICES FOR 
THE SFCOND SYSTEM ~ POINTWISE — AS (3X3) . 


10 


READ 

CONTINUE 


(2) (PH J,K) , K = 1*3) 
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DGMPRD (PI,BT,BS,6,6,3) 


C 

r 

c 

c 

c 

c 

c 


CALL 


*<‘***^^« COMPUTING RESIDUALS ♦♦♦♦♦♦♦♦ # 

*«♦*«♦♦♦<<««««♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦«««<<«♦♦♦♦»♦♦♦ ♦ 

«*«**«** « 


15 


20 

25 


no 15 K 


1 » 6 


KK 

=r 

JJ + K 


A ( KK ) 

= 

O.DO 


KK 

= 

(I-l) ♦ 

3 

00 15 L 

=■ 

1 » 3 


KM 

= 

KM + 1 


A(KK) 

= 

A(KK) ♦ 

BS(K,L) * A(KH) 

DO 20 L 

= 

1 » 3 


LL 


JJ ♦ L 


KK 

= 

LL ♦ 3 


W(L) 

= 

AILL) 


A(LL) 

= 

A(KM) 


A(KM) 

= 

W(L) 


CONTINUE 




9FWIND 


1 


REWIND 


2 


RETURN 




END 
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r . ... 

c ’ , 


c * 

C DAR 04 Y . ♦ 

C . ' ' ♦ 

C ♦ 

r. SL'Br?.niITlNC SETS UP A VARIABLE DIMENSIONED MATRIX IN ♦ 

c p^np^'o SxnPAGF MHOE AS RFCUIRED BY • SSP LIBRARY* ♦ 

r. ♦ 




i. 

c ♦ 


r 


SUr-.RCUTlNE DARRAY(HCPE ,I •J,N,M,S, 0 ) 
implicit 9 FAL*F( A-H* 0 - 7 ) 

DIMENSION sm,om 

M1=N- I 

IF(yrrc_i) ioc,lC0,12O 

icr ij^Mj*^! 

MM=.N^ J + 1 
Dr lie K= 1 *J 

Df^ 110 L = lfl 
IJ=IJ-1 

nr 0 ( \ M ) = 5 ( IJ ) 
on Tf; i^c 
no ij^o 
N'M = 0 

no 13C K=1,J 
DO 1?S L=1,I 
IJ=IJ+1 

NF‘ = NM-*- 1 

i;e suj)=d(nm) 

130 NM = NM-^M 
140 .RFTURN 
FMP 
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c 

r 


r ♦ 

r ♦ 

C MTPY ♦ 

C ♦ 

C ♦ 

f 4;4^4’4c4^4«^t4<4c4c4c«««4c:^«:^3^«4^4c«A«4c4c4c«3^^«««^4c4c«4c4c«4r4c4c»4c«4(*^*44c«4t#4c#4*^«]^:^«*^4c««*^«« 

C ♦ 

C ♦ 

C multiply Twn matrices — FINAL OUTPUT IS A MATRIX, ♦ 

r ♦ 

c ♦ 


C :> 4c 4; :>4i^** + :^***^j>^**4c^#^«****»«**>Ct#4c«* ♦♦♦«♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ ♦♦♦♦«♦♦♦♦ 

C 

c 

c 

SU-^ ROUTINE MTPY( AMT,eMTtMl,M?,M3»CMT) 

IMPLICIT REAL *8 (A«H,0-Z) 

HIM? NS ION AMT(M1,M2) ,BMT(M2»M3) #CMTIM1,M3) 

DO 10 I = 1 , Ml 
• DO TO J = 1 , M3 
CMT( I T J) = 0.0 
nr in L = 1 ♦ M2 

10 rMT(I,J) =r CMT(I,J) + AMTHtL) ♦ BMT(LtJ) 

RCTUt^M 

F^f^ 
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APPENDIX II 


Job Control Cards 
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APPENDIX n 


// (?'^00,].0n) ,CLASS = C 

//STFPl FXPC PRnC = FnRTRANn,PARM='MAP t ID ' , TI ME . CMP = ( 0 , 30 ) 
//CMP.SYSIM no >:>' 


FORTRAN PROGRAM DECK • 


/* 

//STEPP EXEC PRnC=RUMFnRT,PARM.LKED='OVLY, LIST,MAP' , T T M F. LK ED= ( 0 , 20 ) , 

// TIME. On=(3, 10) , REGION. r70 = 252K 

//LKEO.SYSLIB DD DSMAME=S YS 1 . FGRTL I B, D I SP= SHR 

// DD DSNAME=SYS2.F0RTSSP»DISP=SHR 

//LKED.SYSLIN DD DS NAME STEP 1 .CMP. SYSLIN , 01 SP= ( OLD, DELETE) 

// DD « 

OVERLAY ALPHA 


INSERT EULERS, SCALE 

OVERLAY BETA 

INSERT TFnRM,RESIDU,MTPY, SETUP, CSTRNT,DARRAY 


/- 

//GO. F TO IE 001 
// 

//GO. FT02E001 
// 

//GO.FT03F001 

// 

//GO. FT04F001 
// 

//GO. FT07F001 
//GO. SYS IN DD 


DD IINIT = SYSDA,SPACE = ICYL, ( 1, 1) ) ,DISP = ( NEW,DELETE ) , 

DCB=( RFCFM = VBS,l.RECL=600, BLKSI 7.E = 60A) 

DD l!l'IT = SYSDA,SPACE = (CYL, (1,1) ) , D I SP= ( ME W, DEL ET E ) , 

DCB=( RECFM=VBS,LRECL=600, BLKSIZE=604) 

DD t)MIT = SYSDA,SPACE = (CYL, ( 1,1) ) , DI $P = ( NEW, DEL ETE ) , 

DCR=( RECFM=VBS,LRECL=600, BLKSIZE=604) 

DD ll^'IT = SYSDA,SPACE = (CYL,( 1, 1) ) , D I SP= ( NE W, DEL ETE ) , 

DCP=(RECFM=VRS,LRECL=600, BLK size =604) 

DD SYSOIIT = P 


DATA DECK 


/* 

// 
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